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SUMMARY 


This roport pre*icmtn the roaultfl of an analytical and experimental 
denign ntudy aimed at the development of Improved welded metal bellow* 
for application in cryogenic seals. Performance character! ntics for the 
full range of bellows configurationo wore Investigated in the linear and 
non-linear rango of operation. The major characteristic of interest was 
the change in bellows effective diameter over the operating range. Othor 
characteristics of major interest wore the bellows spring rate and the 
maximum operating stresses. 

The analytical work culminated in the selection of opposed and 
nested toroidal bellows plates stiffenod by means of alternating stiffener 
rings as the most promising bellows configurations. Opposed toroids 
appear to have the capability of maintaining constant effective diameter 
within 0. 3 percent for the conditions of interest. Radial corrugations in 
the plates may further improve this. Nested toroids, because more 
convolutions are possible, appear capable of meeting the 0. 1 percent 
goal. Thus, the theoretical work res. ts in the high degree of con- 
fidence that significant advancement may be made in the state of the 
art. 


However, fabrication problems arose concerned with the joining 
of the flexible bellows plates and stiffener rings required for the 
promising configurations. These problems prevented the experimental 
verification of the analytical results within the scope of effort available 
to the program. 

Thus, while the program has not resulted in final experimental 
conclusions regarding performance of the promising designs, it has 
accomplished a significant systematic evaluation of a wide range of 
bellows configurations from the standpoint of their ability to provide 
constant effective diameter. Moreover, it has provided considerable 
experience with the problems invo. zed in joining of mixed bellows 
stacks consisting of flexible plates and stiffener elements where there 
are significant variations in plate thickness among the elements, 
particularly where there are large diameter to span ratios involved. 

This type of structure is felt to be inevitable in achieving substantial 
improvement in the state of the art and, thus, the experience gained 
in this program will be helpful in deve 1 'ping successful fabrication 
techniques. Finally, the test methods „nd apparatus developed and 
built for the accurate measurement of effective diameter under 
cryogenic conditions, and the tooling available for producing experimental 
bellows plates are important tools in the eventual development of 
constant effective diameter bellows. 

It is recommended that these accomplishments be used as the 
basis for a continued program to improve the state of the art in 
bellows seals. Specifically, the next taskB should be aimed at 
developing improved element configurations and joining techniques to 
overcome the fabrication problems. 


INTRODUCTION 


is, *a ffs*£j=MS£sasS3 l r 

the than bo In bollowo oifoetivo diameter over the oporatin iat ^> 0, . 

Other characteristics of major intercot were the bellows npring r<Ao 
and the maximum operating stresses. 

Thn analytical work culminated in the selection of bellows con- 

with StaHw having disparate 

ractto^t£kn“ 8 aa°V prevented experimental verification of the analytical 
results within the scope of the program. 

The orocram did attain significant technical accomplishments 
in term?of P the 8 systematic evaluation of a wide range of bellows con- 
figurations, fabrication experience with the type of structure • 

to 'achieve a conatant affective diameter bellows, and the :J • 
of test methods and apparatus for the accurate meatm* *...t 
effective diameter under cryogenic conditions. 

This work, which it is noped will form the basis for future 
development of improved bellows Beals, is covered in this report. 


Background 

Face contact seals are used to separate the P£?P ell ®“J J T ° m 
the turbine drive gas in rocket engine turbopumps. The useful 
one rating life of these face seals iB highly dependent upon the face 
lnadinc This load must be maintained above a minimum va ^ u ®» 
determine dbv the allowable seal leakage, throughout the operating 
nfessme range of the seal. Thus, the degree of control that can be 
maintained over the face loading throughout this pressure range 
determines the maximum face loading and, consequently, the life of 

the seal. 

The common method for controlling the seal face load is to 
balance the force due to the pressure distribution on the seal face 
with an equal and opposite force due to the pressure acting on a 
secondary metal bellows seal. In theory, since these forces i are bo 
funrtions of the differential pressure across the seal, the seal ia 
load can be maintained constant throughout the pressure range of he 
s'eat, C and is 'determined by the value of the axial spring force present 

in the bellows. 
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This Ideal situation, however, cannot be achieved in practice 
since the effective diameter or area of a bellows can vary considerably 
with differential pressure and operating tomporaturo. Effective diameter 
variations of 1 to 5 percent aro common over the operating pressure 
range within the present state-of-the-art in bellows technology. This 
factor represents a major limitation to the development of improved 
balanced seals for advanced rocket propulsion systems. 

Objective 

The overall objective of the program is to develop configurations 
for welded metal bellows which will minimize changes in effective 
diameter while providing acceptable nominal characteristics throughout 
the desired operating range. The immediate objectives of the study 
reported here are the identification of promising design concepts to 
meet the overall objective, the analysis of these concepts to determine 
their feasibility, and the experimental verification of this analysis using 
special test bellows. 


Method of Approach 

The method of approach used in this study consisted of the 
following major steps: 

Conceptual Design. - The identification and classification of 
promising bellows configurations. 

Screening Analysis. - Linear analysis and preliminary non- 
linear analysis to screen out concepts which are not capable of meeting 
nominal operating requirements or which are not promising for providing 
constant effective diameter. 

Non-linear Analysis. - To predict the non-linear operating 
behavior of the most promising configurations and to select the final 
design configurations for subsequent experimental evaluation. 

This approach provided economy in the analysis of the concepts 
since the most detailed non-linear analysis was performed only for 
those concepts having promise of success. The preliminary analysis, 
however, resulted in a tabulation of the salient features of the whole 
range of concepts and, thus, provides convenient future reference 
material. Furthermore, the systematic approach followed helps 
ensure that no important concept categories were overlooked. 
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description of the system 


A tvnical balanced-pressure, face -contact seal with a welded 
motal hollows secondary seal is illustrated in Figuro 1. The nomenclature 
Tnd terminology for the following technical discussion are presented in 

Table I. 


Considering the axial free -body diagrams of Figure 2 » ***** 
force equilibrium gives the following result for the contact force, 

at the face seal, 


where 


AP 


»■„ • T < D o 2 - D /> A P 

(1) 

P o ■ P i 

(2) 


It should be noted that it is assumed that the nature ° f th ® 10 
pressure distribution across the face seal is known to a sufficien 
degree, either theoretically or by experimental measurement, to 
establish the effective sealing diameter, D For example if a radial 
pressure distribution P (r) exists across tfie face seal, the effective 

sealing diameter, D , can 

8 “ 1/2 


► 

__1 
" A 


can be shown to be given by, 


(P D ‘ 
' o o 


P i D i 


- 8 


r o/ 1 

f P(r) rdr) 

D./2 J 


(3) 


The reaction force, F, , between the bellows and the seal ring 
consists of Tcomponent due K the axial compression of the bellows 
from its free length, and a component due to the net axial pressure 

force acting on the bellows. 

This can be expressed as, 


F b - kA 1 + f (D 0 2 - D.*) AP (4) 

where D is the effective diameter of the bellows. Combining 
Equations (1) and (4) gives the result, 

F s = kA 1 + f- (D, 2 - D e 2 ) AP (5) 

for the contact force on the face seal. 



NAWOHM Utrv. o*fl) 




TABLE I 

Table of Nomenclature and Terminology 
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core 


Effective diameter of the bellows 
Inside diameter of the bellows 
Outside diameter of the bellows 
Effective sealing diameter of the face seal 
Contact force at the face seal 

Reaction force between the bellows and the seal ring 

Bellows spring rate 

Bellows free length 

Axial compression in bellows 

Inside pressure 

Outside pressure 

D - D. 

° i- _ span of the bellows 

Thickness of the bellows plates 

Flexible element of a bellows assembly without 


the end fittings 

convolution = Two bellows plates joined at the inside diameter 
pitch = Axial distance between adjacent weld beads 
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In designing a balanced pressure seal, it is attempted to sot 
D equal to D . Thus, for the ideal case, the contact force will be 
independent of°the differential pressure. The spring force kAf is 
designed to provide the minimum forco necossary for adequate sealing 
withthe minimum operational value of Ai, as determined by the possible 
shaft motion and seal wear. 

In practice, however, even though and D q are balancod with 
no pressure differential across the bellows, the effective d jameter will 
change with pressure differential, A P, unbalancing the seal. It should 
Se noted that! in general, the spring rate, k and the bellows free 
lenffth f will also change as a function of A P. However, if these 
quantities are consideredconstant in Equation (5), then any nonlinearity 
Keen seal force, F„, and 4P can be attributed to a variation in D g 
with AP. This is the 8 convention adopted in this study. 

From a design point -of -view, changes in D have the effect of 
increasing the maximum force which occurs on the seal face, decreasing 
its life. 8 For example, if D increases with A P, then the second term 
of Equation (5) will become negative as A P is increased. This 
requires a larger initial load, kAi , as A P increases. 
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DESIGN GOALS 


The major bollowe design criteria for this program are 
summarized in Table II. These arc discussed as follows. 


The specified ranges for outside diameter, span, and free length 
are commensurate with the space envelope available for the bellows in 
the typical application. 

The m£ -num spring rate specified is consistent with values 
attainable in this envelope within the bellows state-of-the-art. It may 
be seen from Equation (5) that the seal face load is directly proportional 
to the spring rate times the bellows compression, Af, The com- 
pression will vary during operation due to seal wear and axial shaft 
motion. Thus, the seal load variation during operation can be 
minimized by proving minimum spring rate. For a given preload, 
however, the initial bellows compression is inversely proportional to 
spring rate. Thus, the minimum spring rate, while maintaining a 
reasonable value of initial compression, is a desirable design goal. 

The specified axial compression ensures that the seal will not 
become unloaded with expected values of seal wear and shaft motion. 
Pressure capability of 500 psi is sufficient for the turbopump 
application under consideration. The ability of the bellows to withstand 
this pressure differential in both directions would provide maximum 
flexibility in the application of the bellows to various situations. 

However, a bellows design which is capable of providing the desired 
linearity but with pressure capability in only one direction would also 
be useful. 

The 300 cycle life is adequate for the limited life applications 
of interest. 

The seal material specified is Inconel 718 for its excellent pro- 
perties at cryogenic temperatures. 

The ultimate cryogenic fluids will be LH^ and LOg. For test 
purposes, however liquid and gaseous nitrogen will be used. 

The specified allowable variation in bellows effective diameter of 
0, 1 percent ensures that the variation in seal face load due to pressure 
variations will be consistent with the nominal seal load, For example, 
typical values of the spring load, kAi , in Equation (5) are of the order 
of 100 lbs. (444. 8 newtons). 


If D is assumed to change by an amount AD and D is 
initially eqvfal to D , then Equation (5) may be rewritten in me form, 


F = kAf + 


TT 

4 


A P 


A D 

t 

D 


2 

•) + 


AD 


u 


D 


( 6 ) 
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table XL 


Do sign 

Outside Diameter 
Span 


Free Length 


Criteria for the Bellows 

5 to 7 inches 

not greater than 0. 060 inches 
not greater than 1. 5 inches 


Spring 

Room 


Rate (Ap * °)> 

Temperature 


less than 800 lb /in 


Axial Compression 
Pressure Range 
Cyclic Life 

Material 

Medium 

Effective Diameter 


0. 15 inches minimum 

o to 600 psi, both directions 

greater than 300 pressure 
cycles 

Inconel 718 

LH 2 or L0 2 , testing in nitrogen 
Constant within 0. 1 percent 


for the specified conditions of Table II, this results in a possible 
pressure force variation of 40 lbs, (178 newtons). Within the 
existing state-of-the-art, however, a pressure load variation of 400 
lbs. (1780 newtons) is possible resulting in much more serious soal 
face wear. 

It should also be noted parenthetically that the bellows is 
expected to operate over a wide temperature range which alone will 
result in changes up to 0.4 percent in effective diameter and free 
length. Thus, to preserve pressure balancing consistent with the 
above specifications, the ultimate overall seal design must provide 
extremely uniform thermal expansion so that D s and D g change by an 
equal amount. 


BELLOWS PLATE CONFIGURATIONS 

At the outset of this study, a concept generation task was per- 
formed to identify and classify the major bellows convolution shapes 
of interest. Table III shows the range of shapes considered, classified 
according to family and type. 

Within these generic groupings, of course, there are an infinite 
variety of specific configurations and it is not intended to illustrate 
this range. For example, many specific nested ripple configurations have 
been practically applied in the past. Also, useful shapes could be found which 
classify between the types shown. For example, the common single -sweep shape 
probably lies between the nested toroid and the nested ripple. 

It is felt, however, based upon the results of the subsequent 
analytical work that the spectrum of performance characteristics for 
the range of practical shapes is bracketed by these generic types. 

Special mention should be made of the orthotropic* flat plate 
and toroidal convolutions. These shapes resulted from preliminary 
indications that improved performance might result from increasing 
the radial stiffness of the plates without commensurate increases in 
the circumferential stiffness. For the flat plates, this can provide 
increased pressure resistance without corresponding increases in the 
spring rate. For the toroidal configuration, the resulting unbalance 
between radial and circumferential stiffnesses offers the possibility 
of maintaining the toroidal shape more perfectly under axial deflection 
and, thus, a more constant effective area. 

Although the bellows plates were analyzed quantitatively in the 
range of interest as described in subsequent sections, several generalized 
comparisons may be noted here. 


Orthotropic - A specific type of anisotropy where properties 
(in this case, section stiffness) differ along locally 
orthogonal axes. 
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TABLE III 


BELLOWS PLATE CONFIGURATIONS 

Tvdo Convolution Shapo 
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TABLE III (Coat, ) 



Full, Oppoacd 

8 

i 

i 

i 

1 

■ 

fT*^ Mni A *■» 1 

Partial, Opposed 

s 

l 

I 

1 

ioroiaai 

Nested 


1 

i 



Orthotropic (Radial 
Corrugations) 


i 



Full Toroidal 
Segments 

o 

r\^ 

i 

» 

i 

S -Shaped 

Shallow Segments 


i 

! 


Inclined 

a 

i 

• 














A malor distinction may bo drawn between full or nearly full 
toroidal plates and the remaining shapes, in the former the stresses 
duo to pressure are largely membrane stresses and are nearly constant 
thn olnto in the latter, radial bonding stresses are typically the 

Component Id the mognltudo usually varies considerably across 
tho nnnn depending upon the particular shape. The toroid, becauo 
S lts P largely membrane action, has a comparatively high resistance 
to pressure load for typical pinto thicknesses, 

The shallow shapes, however, are typically more compliant 
and more convolutions may be fitted into a given core len B^- Thu0 ' 
they typically exhibit a lower spring rate and a greater 
for\ong stroke. Of these, those plates having a L P flat 

of the span with zero inclination, such as the zero inclination flat 
plate and certain specific nested ripple configurations, have the lowest 

axial spring rate. 

In general, it can be stated that any plate in an °PP ose * °* 
symmetric 8 convolution will exhibit greater linearity (minimum change 
ineffective diameter) than its counterpart in a nested convolution. 

Also, the full toroidal plate is theoretically the most linep" con- 
figuration under pressure load alone. 

The effect of the radial corrugations in the orthotropic plates 
is mainly the reduction in circumferential stiffness. This reduces 
the arfal spring rate of the plates from tne uncorrugated case. Radial 
hendiM stiffness is increased somewhat and the corrugations introduce 
circumVrenttal bending stresses not present in the un-rrugamd case. 
Tt la nossible that, in some situations, the net effect would be an 
increased linearity over the required axial stroke without a reduction 
in pressure resistance. 

The plate and convolution configurations described in this 
section serve as a basis for the detailed analytical work described 
in the following sections of this report. 
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ANALYSIS OF BASIC PERFORMANCE 

Tho design criteria listed In Table II consist of three major 

typo a; 

(a) Dimensional, Material, and Environmental Criteria 

(b) Nominal Peri’ormanee Criteria 

(«) Linearity Criterion 

The dimensional criteria specify the range of bellows diameter, 
span, and length of interest in this application. All of the analysis 
and experimentation performed in this study was concentrated in this 
size range. For the large diameter to span ratios of interest, there 
is no strong variation in performance with diameter. Thus, the size 
of the bellows was fixed at 6 inches (.152 m) mean diameter which 
allows the specified span variation within the specified envelope. In 
order to maintain a low spring rate and provide ease in fabrication, 
it is desirable to keep the span above 0.25 inches. (.0064 m) Thus, 
the extreme cases of 0.25 inches (.0064 m) and 0.5 inches (.127 m) 
were analyzed for their basic performance. The free length specification 
places a limitation on the maximum number of a given convolution 
that can be used, and also, along with the axial compression require- 
ment, sets a minimum of 10 percent axial compression. The plate 
thickness represents the size parameter with the greatest possible 
range and the performance for each type of plate has been studied 
across this range. 

The specified bellows material is Inconel 718 for its excellent 
low -temperature properties. The main effect of the material on the 

analytical results is that it imposes a design stress limit and thus 
limits the range of useful application for a given design. Since the 
linearity criterion is the predominant one, the design stress limit 
has conservatively been set as the minimum Yield Strength in the 
temperature range of interest. The specified 300 cycle cyclic life 
provides, in general, a less stringent stress limitation for Inconel 
718. For reference in the remaining discussion of this report, the 
Yield Strength and Ultimate Strength for Inconel 718 are given in 
Figure 3 as functions of temperature. (1,2)* 

The only environmental condition of interest presently is the 
surrounding medium which has the effect of imposing the temperature 
level at which the bellows must operate. Of the specified media, 

LH 2 has the lowest boiling temperature at atmospheric pressure, 

-423°F (20. 36° F), Thus, the bellows should be capable of operating 
from ambient temperature down to -423 °F. (20. 36° F). 


Numbers in parentheses refer to the list of references. 
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The nominal performance critoria described the capability of 
the bellows to perform their basic function. That is, to provide a 
nominal sealing force and to withstand the combined stresses due to 
axial deflection and pressure differential over a life of 300 cycles. If 
a bellows design is inherently incapable of satisfying these criteria, 
its linearity is of no practical interest. In general, analysis of the 
nonlinear behavior of a bellows is more detailed and costly than linear 
analysis. Thus, the method of approach followed was to screen the 
various types of bellows for their ability to meet the nominal per- 
formance criteria within the specified size range. This screening 
was accomplished by analysis of the linear range of operation, as 
described in the remainder of this section for the various plate shapes 
considered. 

Nonlinear analysis to predict the change in effective bellows 
diameter under load was then performed only for those promising 
bellows shapes, predicted to meet the nominal performance criteria. 


Corrugated Plates 

The general convolution shape for corrugated plates was 
illustrated in Table III. t The opposed plate convolution is of main 
interest in reducing the non-linearity that arises due to differences 
chat arise in the reaction moments at the edges of asymmetric nested 
convolutions under load. For evaluating the capability of a given plate 
shape to meet the nominal performance goals of Table II, however, 
the nested configuration is the most favorable case. This arises 
from the fact that more convolutions can be included within the 
allowable core length, thus, minimizing the spring rate and the stresses 
due to the specified axial compression. For this reason, the nested 
configuration was analyzed in this study. 

Because the nested ripple plate is in common use as a bellows 
seal element, considerable performance information has been 
accumulated. For example, it is known that the predominant stress 
components in these elements are the radial bending stress and the 
tangential direct or membrane stress. The distribution of these 
stresses is quite dependent upon the diameter to span ratio. For 
small values of this ratio (i. e. up to about 10), the magnitudes of 
these components are comparable and the distribution follows the 
plate ripples, often having large amplitude of variation across the span. 
In this range, considerable variations in performance can occur as 
the detailed shape of the plate is varied. 

For large values of the diameter to span ratio, which is the 
range of interest in this study, however, the bending stress becomes 
predominant as shown in Figure 4. This approaches a straightline 











variation across the span which is characteristic of flat plate 
behavior. In this range, the major factors that affect performance 
are the overall plate inclination and the degree of coning or 
inclination near the edges rather than the details of the ripples near 
the center of the span. 

For this reason, the analysis discussed in this section is 
based upon a standard plate shape developed by Sealol Inc. for seal 
applications. It is felt that this represents a typical well -designed 
bellows plate of the nested ripple configuration. 

The solution of Figure 4 and those upon which the remaining 
data of this section is based were obtained using the computer program 
described in Appendix A. This computer program solves the generalized 
shell equations as applied to bellows elements in the linear range. 

It, thus, provides an approximation to the nominal performance of 
a bellows. Significant errors would be expected for the actual mag- 
nitudes of the stresses under loads which produce large deflections 
and, of course, this program is not capable of predicting nonlinearity 
of the bellows. However, the predicted spring rates have been found 
empirically to hold over considerably large deflections. 

The nominal performance of the tested ripple plates in the 
size range of interest is shown in Figures 5 and 6. In Figure 5, 
a mean diameter of 6 inches (. 152 m) and a span of 0.25 inches 
(. 0064 m) is assumed. The curves plot the peak maximum stress 
intensity as a function of plate thickness. This stress intensity 
combines the effects of a 500 psi (3.445 x 10° newtons/W 5 ) pressure 
load and an axial stroke of 0.15 inches (.00381 m). Figure 6 plots 
the same thing for a span of 0.5 inches (.0127 m). 

Two curves are shown for the nested ripple plate in Figure 5. 

One represents a bellows core consisting of 30 convolutions regardless 
of plate thickness. This would provide sufficient stroke within the 
allowable length of the core only for plate thicknesses up to 
approximately 20 mils (5.08 x 10 ~ 4 m). The other curve allows the 
number of convolutions to vary according to the maximum that can 
be fit into the core length for each plate thickness between 10 and 
20 mils (2.54 x 10-4 an d 5.08 x 10" 4 m). Below 10 mils (2.54 x 10 m) 

the number was held at 60 which is felt to be a practical, although 
arbitrary limit. 

In both figures, the performance of the corresponding flat 
disk bellows is shown for comparison. This will be discussed in the 
next section. 

Figure 7 shows the bellows spring rates for each of the cases 
given in Figures 5 and 6 as a function of plate thickness. The 

discontinuity in the 60 convolution case of Figure. 7 results from the 
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fact that, with a plate thickness beyond 10 mils (2,54 x 10” m) 60 
convolutions can no longer be fitted within the allowable core length, 

Thus, beyond this point, the maximum number which can bo fitted 
in is used. This number, of course, changes as thickness is 
increased. 

The limitations in the capability of nested ripple plates for 
meeting the nominal performance goals can be seen by considering 
the data shown in Figures 5, 6 and 7. Remembering that the 
maximum allowable stress intensity is approximately 150,000 psi 
(1,35 x 10* newtons/m 2 ) as shown in Figure 3, it can be seen from 
Figure 5 that the plate thickness for a 0. 25 inch (. 0064 m) span 
would have to be approximately 9 mils (2. 285 x 10 " 4 m) or larger. 

Figure 6 shows that the 0. 5 inch (. 127 m) span is virtually precluded 
from this application since the stress does not get down to the allowable 
range until the plate thickness reaches 20 mils (5,08 x 10"’ m). No 
further reduction in stress can be achieved through additional con- 
volutions since 30 convolutions is the maximum that can be fit into 
the available space for this thickness. 

Figure 7 and Figure 5 show that, for the 0. 25 inch (. 0064 m) 
span, a spring rate of 800 lbs/in (1.4 x 10^ newtons/m) or less can 
be obtained in the acceptable stress range only for plate thicknesses 
between 9 and 11,5 mils (2. 285 x 10 "4 m and 2.92 x 10 "4 m). 

Thus, it may be concluded that the nominal performance goals 
could be met only for the minimum practical span, 0. 25 inches 
(.0064 m), with a plate thickness of approximately 10 mils (2.54 x 
10 m) and greater than 30 convolutions. This is a very narrow 
range of acceptable parameters and consequently must be concluded 
to be a marginal situation. Any slight variation from the predicted 
values of either stress level or spring rate could eliminate the acceptable 
range entirely. For example a 27 percent error on the low side in 
the predicted stress value, which is easily possible for the methods 
used, could preclude the use of the 0.25 inch (.0064) span bellows. 

For this reason, it is not felt that the nested ripple bellows 
configuration is suitable for this application. 


Flat Disks 

As noted in the preceding section, membrane stress effects are 
small for the corrugated plate shape in the diameter to span ratio 
range of interest. This implies that the corrugations and plate tilt 
have a relatively low importance in this range. For this reason, 
flat disks were analyzed for their basic performance. 
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Accurate relationships are available in the literature tor' the 
linear behavior of flat disks (3). Solutions are also available for the 
behavior of coned flat plates (i. o. Belleville washers) ( 4 )^ ^he® 
are somowhat unwieldy however, and it was found that, for me t rg 
diameter to span ratios of interest, the results aro a WWjdn»ted 
within 10 porcent by considering the plates to act as a wide g 

beam with the same span, thickness, and edge conditions. This is 
sufficient accuracy for the purposes of this screening analy . , 

these approximate^ formulae^ which are summarized below were 

used to obtain the flat disk data plotted in Figures 5, 6 and 7 
comparison with the nested ripple plates* 

A schematic diagram showing the ini tial p re c omp r c s Sion and 
pressure loading for the straight beam model of the flat disk bell 
is given in Figure 8. 


The force per unit circumference, 
axial deflection, A L, is given by, 

F = K A L 


F, required to produce the 


(7) 


where the spring rate per unit circumference, K, may be approximated 

by, 

l3 


K = 


Et" 

S 3 (1 - V Z ) 


(8) 


where 


E 

v 

t 

S 


Young* s Modulus 
Poisson* s Ratio 
plate thickness 
span 


The maximum bending stres. produced by this initial precompre.sion 
may be approximated by the relation, 


a f “ 


_ 3 F S 

.2 


If the bellows is then pressurized with pressure, P, an additional 
bending stress is produced according to the relation, 


(9) 
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(10) 


°p R 


p gf 

2 t* 


Superposition of the two stress components given by Equations (9) 
and (10) at the point of maximum stress whore they add was used to 
calculate the behavior shown in Figures 5 and 6. The spring rates 
shown in Figure 7 were calculated using Equation (8). 

. The salient conclusions that may be drawn from these results 
are that the flat plate bellows have a lower spring rate for the same 
number of convolutions than the corrugated plate hollows. Also, the 
bending stress due to axial compression alone is lower. However, 
the total stress due to both axial compression and pressure is > 
significantly higher than that for the corrugated plates. For P late 
thicknesses where the stress due to pressure is acceptable, the spring 
rate is increased to the point where no advantage is gained in com- 
parison witl the corrugated plates. Thus, the flat plate bellows is 
also concluded to be unsuitable for this application. 

Finally it should be noted that the performance superiority 
shown in Figures 5 and 6 for the corrugated bellows does not con- 
tradict the observation, made in the previous section, that corrugations 
near the center of the span have little effect in this size range. 

Rather the advantages of the corrugated plate are believed to be due 
to the more favorable edge conditions produced by the corrugated con- 
figuration and the resulting smaller effective span. 


Orthotropic Plates 


As depicted in Table HI, orthotropic plates contain corrugations 
around their circumference, lined up along the plate radii* 

As noted in the preceding section, the major drawback for flat 
disks was that increasing the plate thickness to obtain pressure 
resistance resulted in spring rates which are too high. 


This can be mitigated by radially stiffening the center of the 
span while leaving the region near the edges of the span, m08t 

of the bending under axial deflection occurs, unstiffened The 
stiffening of the central part of the span reduces the reac ^°“ 
and the resulting stress due to pressure loading as well ^ as the , 8k J e88 
in the central portion of the span. At the same time, the bending _ 
stress and spring rate under axial deflection are only slightly increased 
since the central part of the span primarily tilts with little bending 
under axial deflection. 


Radial corrugations are a convenient way to achieve this 
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stiffening using stamped plates, and they allow for the possibility of 
nosting the plates in a stack. The radial corrugations also provide 
a decreased circumferential stiffness which aids in the desired action. 

As for flat plates, this configuration was analysed using the 
straight beam approximation, The schematic diagram of Figure 9 
shows tho model used. The angularly flxod boundary condition is 
imposod as for tho flat plato by tho convolution symmetry. 

It is assumed that the center of tho span is corrugated to 
provide a section moment of inertia, I, which is n times that of the 
uncorrugated plate, I . It can be shown from the properties of the 
sections that, ° 


where 

h * amplitude of the corrugations 

t = plate thickness 


Thus, a moment -of -inertia ratio of 100 can be achieved for a 
corrugation amplitude to thickness ratio of four. 


Formulae for the spring rate, K, and the stress due to pressure, 
<Xp, may be computed from energy methods using the method of 
Gastigiiano. The relation for spring rate per unit circumference is, 


K 


Et“ 


S 3 (1 - V 2 ) 



1 

-12 x 2 + 6 x) (1 



( 12 ) 


which is similar in form to that of Equation (8) for fiat plates with 
the addition of the correction term. In this relation, x represents 
the fraction of the span which is uncorrugated at each edge of the 
span. 


Tt*e stress due to the axial compression remains as given by 
Equatio?, ,9) of the preceding section where the force, F, is computed 
from Equation (7) using the spring rate given by Equation (12). 


by, 


The maximum stress due to_a pressure loading, P, is given 

,2 


- P S‘ 
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which again is of the form of the flat plate relation, Equation (10), 
with a correction term. Fur the typically large valuea of n of Internet 
for this configuration, Equation (13) roducee to, 

<r S (3 X - 2 x 2 ) (14) 

P 2 t. 

Equations (7), (9), (12) and (14) wcro used to predict the performance 
shown in Figures 10 and 11. The number of convolutions assumed 
for these solutions was 30. Since there is no experience in stacking 
this shape, this number is open for argument. However, the 
minimum number used for the solutions of the previous shapes was 
30 and this number represents the maximum number that can be 
fitted in with a plate thickness of 20 mils (5.08 x 10 “4 m ). It is 
probably a reasonable maximum, even at the lower plate thicknesses, 
due to the fact that stiffening rings would amost certainly need to be 
provided with this shape. 

Comparing the stresses of Figure 10 with the corresponding 
flat plate solutions of Figure 5 and 6, it can be seen that the radial 
corrugations indeed improve performance in this respect. For example, 
a 50 percent reduction in stress results for the 0*25 inch (.0064 m) 
span with a plate thickness of 10 mils (2.54 x 10 “4 m). For a 
0.5 inch (.0127 m) span, the reduction is even more pronounced, 
almost an order of magnitude. However, the reduction in stress 
is achieved only at the expense of a considerably increased spring 
rate as may be seen from a comparison of Figures 7 and 11. The 
net result is that, based on the nominal characteristics of Figures 
10 and 11, there is a small useable range of plate thickness for 
each span. This is approximately 2 mils wide for the 0. 25 inch 

(.0064 m) span and 5 mils wide for the 0.5 inch (.0127 m) span. 

Above this range, the spring rate is above the 800 lbs/4n (1.4 x 10 
newtons/m) limit of Table II; below it, the stresses exceed the 
150,000 psi (1.35 x 10** newtons/m 2 ) limit. 

Although the nominal stresses are indeed lower than those for 
both the nested ripple and flat disk plates of Figure 5 in this thickness 
range, the design is still felt to be extremely marginal. In no case 
does the predicted stress level fall below approximately 80,000 psi 
(0,72 x 10^ newtons /m 2 ). Local stress concentrations near edge 

fixtures and at the transition points between the radial corrugations 
and the unstiffened edge region could easily raise the stress beyond 
the 150,000 psi (1.35 x 10 ^newton s/m 2 ) limit. 

Finally, the reduced circumferential stiffness of the plates 
due to the corrugations will require stiffening rings to achieve the 
desired degree of linearity just to resist diametral growth under 
pressure. Although this is a common requirement for most of the 
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Figure 10. Combined Stress for Orthotr opic Plate Bellows, 
(Radial Corrugations) Due to Pressur e and Axial Deflectio n 


















shapes considered, it becomes moat critical from a fabrication point 
of view when a large number of convolutions are required as in this 


case. 


Thus, the net cffoct of the orthotropic flat plate bellows is to 
provide a marginal design similar to that of tho normally corrugated 
plate only with a smaller plate thickness. For this reason, it was 
not carried further for evaluation in the nonlinear range, 


Full Toroidal Plates 

The use of toroidal convolutions is well-known for high pressure- 
carrying capability. The welded version would consist of two opposed 
semi-toriods, welded together at one edge and connected to adjacent 
convolutions at the other edge by stiff spacer rings. 

Solutions for the behavior of toroidal convolutions are given in 
the literature. (3,5). The schematic diagram of Figure 12 shows the 
major variables of importance. 

The following formulae may be used to compute the deflection 
and stresses under both precompression and pressure loading in a 
manner entirely parallel to that discussed in the preceding 
Sr Tat disks and orthotropic plates. The formulae are limited to 
regions where the following inequalites are satisfied among the bellows 
dimensional parameters (5). 


4 - > 40 


(15) 


2 t R 


( 16 ) 


In this region, the spring rate per unit circumference is given by (6), 


K = 



10.88 R S 



(17) 


The maximum bending stress due to a force per unit circumference, 
F, is, 
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The maximum stress due to pressure loading, P, becomes, 


°P ~ 


= 2JL 

Zt 


The latter stress is a membrane stress and the approximate result 
applies within 5 percent in the size to span region of interest in this 

study. 

Although stability and buckling of be 'ws elements is always 
a consideration in bellows design, it is of s^ cial concern for toroidal 
elements which often require only relatively small plate thicknesses 
for pressure resistance. When externally loaded, these elements 
could undergo large deformations due to buckling. Since the Present 
application has as a goal the capability of pressure loading in both 
direction, buckling limits were studied to determine whether they 
impose practical design criteria. 

No previous analysis was found for the local buckling of a 
toroidal convolution. However, it is obvious that circumferential 
effects and complex curvature will make the toroid more stable than 
its corresponding straight tube. Available solutions for a straight 
tube with external pressure can be used to set an upper limit on the 
required plate thickness, t, as follows (3). 

(-r> 3 > ^ 

of 30°“ ft P.M2. S? 1 * 6 

10 11 newtons/m 2 ) into this relation gives, 


— > 0. 04 
r 


A more conservative mode of collapse is the local buckling of the 
crest of the toroid. Using, as a conservative approximation, the 
buckling of a 90° straight arch, a lower limit for the required plate 
thickness (3), 


(J_\ > = 

\ * TT! 


0.02 
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for the above numerical values. Thus, the range for realistic buckling 
limits lies between the values of Equations (21) and (22), Above this 
range, tho toroid is almost certainly stable and below, it is likely 
to be unstable, 

Tho relationships of Equations (17) through (19) were used to 
prodict the combined stresses for toroidal bellows having spans of 
0,25 inches (.0064 m) and 0,50 inches (.0127 m) in Figures 13 and 14, 
respectively. The stability limits discussed above are also indicated 
on the graphs, The spring rates for these bellows are shown in 
Figure 15. 

It can be seen from Figures 13 and 14 that the combined stresses 
for both spans are within the Yield Stress for the material (150,000 psi 
or 1. 35 x 10 4 newtons/m 2 ) in the thickness range between 3 and 10 
mils (7.62 x 10" 5 m and 2.54 x 10 -4 m). The stresses are most 
favorable for the 0.25 inch (.0064 m) span where more convolutions 
can be accommodated within the core length and the pressure stresses 
are smallest. This span, of -course, also provides less restrictive 
stability limits. Figure 6 shows that virtually all combinations of span 
and number of convolutions give a spring rate within the design goals 
in the plate thickness range of interest. 

The location of the peak combined stress is also of interest for 
the toroid. The stress due to pressure loading is nearly constant 
throughout the toroid. The bending stress due to axial loading is 
maximum at an angle of approximately 45® along the toroid from its 
mid span. Thus, the maximum combined stress does not occur at the 
edges of the span where stress concentrations due to the weld joints 
may be expected. 

Modification of the toroidal plates to add radial ribs similar 
to those studied for the orthotropic flat plates was investigated. The 
purpose of this modification in the toroidal plates was the mitigation 
of shape changes which tend to bring about nonlinearity under axial 
deflection, and, as such, will be discussed in a subsequent section 
treating nonlinear considerations. It should be mentioned here, however, 
that in the orthotropic toroidal bellows, the nonuniform stresses caused 
by axial deflection of the bellows are expected to be relieved in addition 
to the provision of greater linearity. 

In summary, the full opposed toroid appeared to offer basic 
performance which was within the desired design goals. From a 
practical design point of view, fewer convolutions were required and 
did not appear to present outstanding fabrication difficulties. Although 
stiffening rings would be required to achieve the deBired degree of 
linearity, their number and cost is minimized due to the small number 
of convolutions. Thus, it was recommended that this configuration be 
further evaluated as one of the most promising for this application. 
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Figure 15. Spring Rates for Toroidal Bellows 


Partial Toroids 


The advantages of tho toroidal shell for pros sure -carrying 
capability are apparent from tho results of the provious section. These 
advantages are due to the membrane action which only depends upon 
the fact the radius of curvature of the shell is constant across 
the span. Except for boundary effects near the points of attachment 

similar advantages would be expected for partial toroids where each 
plate consists of a section including less than 180°. These have 
the advantage that more convolutions could be fitted into a given core 
length if the plates were nested rather than stacked symmetrically. 

Thus, the basic performance of the partial toroid was investigated 
as described in this section. 

Partial toroids can be assembled in various configurations, both 
with and without stiffening spacer rings as shown in Figure 16. The 
angular stiffness of the spacer rings is extremely important, 
especially for the nested configurations. 

' Qualitatively, the performance of the partial toroids can be 
understood by considering .the normal toroidal membrane action coupled 
with the resultant reaction forces produced at the boundaries or edges 
of the span. The membrane stresses due to pressure loading act 
along a direction tangent to the plate. Thus, in Figure 16(a), the 
reaction forces, Fa and F B , are almost equal in magnitude but have 
a small resultant force, Fp, due to their difference in direction. 

This reaction force tends to tilt the convolution and expand it radially 
away from the centerline. The effects of the reaction force which is 
always present for the configuration where the plates are joined 
directly to each other can be reduced by providing a stiff spacer 
ring to resist the deformation as in Figures 16(b) and (c), This will, 
of course, produce local bending stresses near the stiffening rings. In 
Figure 16(d), the convolution is oriented so that the tilting moments 
due to the reaction force at each end cancel each other. The con- 
volution, however, would still tend to deflect radially outward. 

By attaching the plates to either side of the stiffening ring as 
in Figure 16(e), the plate tangents, and thus the direction of the 
resultant forces can be made parallel. In this way the reaction 
loading on the stiffener can be reduced to a moment or couple which 
would have to be resisted by the stiffener if tilting is to be avoided. 

If the plates are stacked in a symmetrically opposed configuration 
as shown in Figures 16(g) and (h), there can be no section tilting but 
the resultant force, F«, tending to deflect the convolution radially 
is much larger. This can be resisted, of course, if sufficient 
radial stiffness is provided by stiffening rings as in Figure 16(h). 
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Figure 16. Partial Toroid Convolutions Having Var ious Configuration:; 


Methods of Assembly, 










The partial toroids wore analysed In detail, using tho linear 
computer solution described In Appendix A, for two limiting configurational 

(a) The caao where the plates are attached together by 
moans of sufficiently stiff spacer plates so that the 
span edges can be considered angularly and radially 
fixed, 

(b) Tho simple nested configuration of Figure 16(a) with 
no auxiliary stiffening. 

The variation of the total combined stress and spring rate for 
six convolutions as a function of the arc semi -angle, ^ , are plotted 
in Figure 17 for the fixed edge conditions uith a plate thickness of 
5 mils (1,27 x 10 m ). The results show that the combined stresses 
are acceptable for this whole range of angle, ^ , with a minimum 
in the range of 70°, 

The variation of performance with plate thickness for a given 
arc segment ( $ = 60°) is shown in Figure 18. Here it can be seen 
that the stresses are at a minimum in the range of 3 to 4 mils 
(7.62 to 10.16 x 10 "5 m ) plate thickness. 

Both Figures 17 and 18 show that the maximum stress 
intensity and the spring rate for six convolutions are acceptable 
within the design goals for this program provided stiffening rings are 
used so that the edge condition for the bellows plates approach the 
angularly and radially fixed case. The plates could be stacked 
symmetrically opposed as shown in Figure 16(h) or nested as shown 
in Figures 16(c), (e) or (f). In the latter case, many more than 6 
convolutions could be used, further reducing the spring rate and the 
stress due to axial stroke. 

Nested partial toroids without stiffening rings as shown in 
Figures 16(a) and (d) were also analyzed. The results showed that 
severe shape distortions occur under pressure loading and the combined 
stresses are intolerable. 

Referring to Figures 16(a) and (d) for the definition of 


variables, 

cases 

were 

run with: 


-e- 

ii 

60° 



a = 

0. 1 

radians 

o < 

P < 

0. 1 

radians 


41 


Intensity, (psi x 10 (newtons/m x 1. 45 x 10 


Combined Stress due to Pressure and 

Axial Deflection 

1 / 


N ^Spring Rate 


, Mean Diameter - 6 in (. 152 m) 

Span - 0. 25 in (. 0065 m) 

Core Length - 1. 5 in (. 038 m) 4 

Plate Thickness - 5 mils (1. 27 x 10 pa) 
Number of Convolutions - 6 6 

Differential Pressure - 500 psi (345 x 10 

newtons /m ) 

Stroke - 0, 15 in (. 0038 m) 


Arc Semi -Angle, 0 , Degrees 

Figure 17. Variation in the Performance of Partial Toroids with Arc Lejigth 

For Fixed Edge Conditions 


Spring Rate, lbs/ins (newtons /me ter x 0. 571 






Stress Intensity, (psi x 10 ) (newtons/m 2 x 1. 45 x 10 


Spring 

Rato 


Combined Stress due to 
Pressure and Axial 
Deflection 


Mean Diameter 6 in (. 152 m) 
Span - 0, 25 in (. 0065 m) 

Core Length - 1. 5 in (. 038 m) 
Arc Semi -Angle - 60 degrees 
Number of Convolutions - 6 
Differential Pressure - 500 psi 
(345 x 10° newtons/m^ 
Stroke - 0, 15 in (, 0038 m) 


Plate Thickness, mils (in x3, 939 x 10 ) 


Fieure 18, Variation in the Performance of Partial ToroidB with 


Plate Thickness for Fixed Edge Conditions 


Spring Rate, Ibs/in. (newtons /meter x 0. 57 x 10 






A typical convolution distortion is shown in Figure 19. Variation of 

tho mean cone angle, , was found to have little effect on the 

stresses. The major stresses, on tho other hand, are approximately 
in direct proportion to the pitch angle. 

For a span of p,. 25 inch (. 0064 m) and a pitch angle, , 
of 9.1 degrees (pitch -50 mils (1,27 x 10-3 m)), the combinod stress 
intensity predicted is 400,000 psi, (2.76 x 10? ncwtons/m*) with a 
pressure differential of. 500 psi (3,445 x 10o newtons/m^). This 

assumes 10 convolutions and a stroke of 0.15 inches (3.81 x 

10-3 m ). 

In order to reduce the stresses to the level of those for the 
stiffened ring as shown in Figure 18, the pitch would have to be 
reduced by a factor of about 5, That is, approximately 50 convolutions 
would be required. For this case, plate clearance would be a problem 
and fabrication would be costly. Although span and thickness could 
be increased to reduce the stress or allow fewer convolutions, in- 
dications are that the basic performance would not be significantly 
better than the corrugated or flat plate shapes. 

In summary, the partial toroid plate shape are attractive if 
stiffener rings are used to provide angular and radial fixity at the 
edges. This configuration is recommended as the most promising 
nested design. Without the stiffener rings, the partial toroid does 
not appear attractive for this application. 


S-Shaped Plates 

The final family of plate shapes shown in Table III to be 
analyzed for basic performance was the S-Shaped family. Preliminary 
consideration indicated the possibility that, in the nested con- 
figuration, the reaction moments of successive plates would tend to 
cancel due to the difference between the sign of the curvature for the 
two plates with respect to the pressure load. This would provide a 
convolution resistant to the large convolution deformation of the type 
shown in Figure 19. When stacked symmetrically opposed, it was 
not expected to offer any advantages over the previous configurations 
analyzed. 

The dimensional parameters used to define the S-Shaped plates 
are defined in Figure 20. The range of these parameters considered 
were as follows: 
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Figure 20. Definition of Dimensional Parameters for the 

S- Shaped Plates, 
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Plate 


-4 

Thickness, tj S mils and 7 mils (1*27 x 10 m and 
1. 778 x 10-4 m ) 

Cone Angle, a t 90°, 120°, 180° 

Wave Length, \ : -jr t° -sjjr y2 

Arc Angle, <f> : 90°, 120°, 180° 

Span, S: 0.25 inches (.0064 m) 

Mean Diameter, Dj 6 inches (. 152 m) 

The various shape classes formed by varying these parameters 
are shown schematically in Figure 21. 

These shapes were analyzed for stress and deflection under 
both a 500 psi (3.445 x 10° newtons/m 2 ) pressure loading and an axial 
stroke of 5 mils (1.27 x 10 " 4 m) per convolution. This stroke is 
consistent with achieving the specified bellows stroke of 0.15 inches 
(3.81 x 10-3 m ) if 30 convolutions are used. This, of course would 
not be possible unless the plates were nested but represents the most 
optimistic assumption. 

The stresses and spring rates were computed with computer 
program, NONLIN, developed by Battelle Memorial Institute (6). This 
computer program which is described in Appendix 3 is capable of 
predicting the nonlinear behavior of a bellows and was also used in 
this study for that purpose as discussed in a later section. It is also 
capable of predicting the linear performance, however, and is con- 
veniently designed for ease in treating various plate shape inputs. 

Thus, it was used to investigate the basic performance of the 
S -shaped plates. 

The plates were assumed to be angularly fixed at the edges. 

This is again an optimistic assumption. It represents closely the case 
where the plates are stacked symmetrically opposed but would require 
stiffening rings for a nested configuration. However, it will be seen 
that, even under the optimistic assumptions, the pressure stresses 
for these plates are too high. Thus, their analysis was not carried 
further under more realistic assumptions. 

The results of the computer analysis are plotted in Figures 22 
and 23 where the results for the nested ripple and the flat disk shapes 
are also given for comparison. Figure 22 shows the predominant 
stresses as a function of plate thickness for each of the plate shapes 
treated under two separate loadings: 
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Ficure 22. Stresses in S-Shaped Plates Compared with Corrugated and 


Flat Plate Shapes 
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/ 2 

(a) 500 psi (3.445 x 10 newtons/m ) pressure loading with zero 

deflection, 

-4 

(b) 5 mil (1. 27 x 10 m) per convolution deflection with zero 

pressure load. 

It can be seen that the general behavior is similar to the nested 
ripple and flat disk bellows. The stresses for the S-shaped plates under 
pressure are higher than those of the nested ripple plates and fall on 
either side of the flat disk shape depending upon detailed configuration. 

The stresses due to deflection are all greater than those of the nested 
ripple plate 3 and flat disks. 

The spring rates, given in Figure 23, also have a similar 
variation with thicknesses to the nested ripple plates and flat disks, but 
are much higher, even assuming the optimistic 30 convolutions. 

It can easily be seen that, for plate thicknesses where the spring 
rate approaches the acceptable 800 lb/in (1.4 x 10^ newtons/m) range, 
the stresses in the S-shaped plates are far above the design stress limit 
of 150,000 psi (1.35 x 10? newtons/m 2 ). Thus, the basic performance 
of the S-shaped plate does not warrant its further consideration for this 
application. 

It might be noted parenthetically that Figures 22 and 23 provide 
a convenient reference, summarizing and comparing the basic linear 
performance of bellows plate shapes where bending stresses predominate 
over membrane stresses in the large size to span range. Although 
these families of plate shapes are not promising for the present high 
pressure application, they are generally useful as bellows seal elements 
in less demanding situations. 

Summary of Basic Performance 

The basic performance of bellows plates includes their pressure 
carrying capability and force -deflection characteristics. The analysis 
performed in this study shows that there are basically two categories of 
such performance characteristics: 

(a) Bending Plates - where the plates are relatively flat across 
the span including multiple waves about the chord line. 

(b) Membrane Plates - where the plates form a significant 
part of a semi -toroid. 

The predominant stress in the bending plates under pressure is 
the radial bending stress. For a given span, this stress is relatively 
large and the plates must be thick to withstand it. This, in turn, leads 
to large spring rate and deflection stresses or, conversely, a large 
number of convolutions. 
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In the diameter to apan rango of interest, diameter was not 
found to ha vo a strong effect on performance. 

Both categories were found to be susceptible to large angular 
deformation or coning of the convolution section under pressure loa i g 
when The td ate s are nested rather than being stacked symmetrican y 
onnosed This is most serious for the thinner membrane plates and 
for cases where the pitch is large and requires, in general, the use 
of stiff spacer rings to withstand the load. 

The major conclusion drawn on the basis of the analysis of 
basic performance is that the membrane category offers the bes«. 
combined ability to meet the pressure and spring rate requirements 
of Table II within the specified dimensional constraints. Th ® 86 
eitlwsr be symmetrically opposed full toroidal plates or rested partial 
toroids with stiffening rings. 

The corrugated plate bellows was the best of the bending plate 
category It appears capable of providing marginally acceptable 
performance with a relatively large plate thickness and a large number 
of convolutions. However, its performance ts inferior to both of the 
promising toroidal configurations and, thus, it was not evaluated 
further for this application. 
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NONUNEAR CONSIDERATIONS 

Tho moat promising bollows plate configurations for mooting 
the pressure and spring rate requirements of Table II within the 
specified dimensional constraints wore found to be symmetrically 
opposed toroidal plates and nested toroidal platos, each with stiffening 
rings. In the following sections, their linearity or stability of 
effective diameter is studied analytically in order to select tho most 
promising test configurations. 

The design goal for the bellows is that it will maintain its 
effective diameter constant within 0. 1 percent at a given axial com- 
pression or stroke over the entire pressure range, 0 to 500 psi 
(3.445 x 10^ newton/m^). 

The method of approach followed was to conduct parametric 
studies to predict the effects of the various design parameters for 
each shape on the change in effective diameter. The optimum com- 
bination of parameters was then selected for further consideration in 
fabrication and experimental studies. 


Symmetrically Opposed Toroidal Plates 


The basic element of the symmetrically opposed toroidal bellows 
is shown schematically in Figure 24 where the principal nomenclature 
for this study is defined. Because of the symmetry, the plate may be 
assumed to be angularly rigid or fixed at the edges as shown. Also, 
it was determined in the preliminary study that sufficiently rigid 
stiffener rings can be provided so that the edges may also be con- 
sidered to be radially fixed. Thus, the model used for parametric 
analysis as illustrated in Figure 24 is free to deflect axially under 
an applied loading but is fixed radially and angularly at the edges. 


In addition to the parameters defined in Figure 24, the following 
variables will be used in the discussion. 


(a) 


N p . by 


Effective Diameter 

Defining an axial bellows force per unit circumference, 






(23) 
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t 



p - pressure 

5 - axial deflection 

S - span 

t - plate thickness 

r - radius of curvature 

- inner radius 

r - mean radius 

m 

R - outer radius 

* ss included angle of the toroidal segment 


Fiflure 24. Schematic Diagram of the Symmetrical ly Opposed 
Toroidal Plate Element and the Definition of Principal Nomenclature. 
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where 


j? n reaction force between the bellows and the seal ring 
b 

D B effective sealing diameter of the face seal 

8 

and combining with Equation (4) allows the definition of an effective 
diameter aB follows, 


D 2 = D 2 + (24) 

e “ o P P 

where P is used to denote the differential pressure across the bellows. 

(b) Nonlinearity of Effective Diameter 

A small nonlinear perturbation, A D , of effective 

diameter about its nominal value, D , may be related to a perturbation 

in bellows force, AN , by differentiating Equation (24). Assuming 

that KAi and D are** constant, this gives, 
o 

2D AD = - p — A N (25) 

e e P p 

which may be used to define a percentage change or nonlinearity in 
effective diameter. Defining the percentage change in effective 
diameter as, 




| A D 

I e 

D 

1 e 


x 100, 


and combining with Equation (25) gives, 


where 


P 

D 


% A D e = 


2 D 


P D 


AN x 100 
P 


A N 

E 

15 


a 500 psi (3. 445 x 10^ newtons/m^) 

a D = 6 inches (0, 152 meter) 

8 


(26) 


It (should be rr.*ed tftat the above definition aflfiumefi that (spring force, 
K4i, la a con a taut under pronfsuro load, In other worda, any change 
that might occur In the (spring load under preaoura in lumped into 
4N and thua, by definition, into AD , This in in accordance with 
the ‘convention dincunnod under "DoncHptson of the byntom, " Thin 
convention in roaoonable from a donlgn point-of -view becauno any 
change in spring force with pressure appears an an effective nonlinearity 
between pressure force and pressure, or an a change in effective 
diameter. Moreover, thero is no method for separating the various 
contributing effects to a nonlinearity detected in a solution for N 
at a given pressure level and initial deflection. p 

Thus, Equation (26) was used to interpret the nonlinear computer 
solutions in the following discussion. 

The only parameter fixed for the parametric study of non- 
linearity was the nominal or mean diameter, D , which was taken 
to be 6 inches (0, 152 meters). The other pararneters were varied 
in the ranges shown in Table IV. 

A total of 50 computer runs using computer program NON LIN 
were made to study these parametric variations as well as the effects 
of the boundary stiffness and certain shape modifications. The 
detailed numerical results from these computer runs are tabulated 
in Appendix D. Typical stress distributions are also presented in 
Appendix D for plates having 60°, 120°, and 180° included angles. 

The effects of the variation of each of the major parameters on change 
in effective diameter over the desired pressure range are discussed 
in the following paragraphs. 

Effect of A xial Deflection. - The percentage change produced 
in the effective diameter between small pressure loads and a pressure 
loading of 500 psi (3.445 x 10° newtons/m 2 ) as a function of the axial 
deflection of a single plate, 6, is shown in Figure 25 for various 
included angles. (Reference Run Nos. 1, 2, 3, 4, 15, 18. 19, 20. 

21, 22 and 34 of Appendix D). 

It can be seen that the percentage change in effective diameter 
is approximately linear with the axial deflection for a given included 
angle and that there is a marked change from a large included angle 
approaching a full semi -toroidal and a small toroidal section. For 
the large included angles, it can be seen that to achieve the desired 
goal of 0. 1 percent or less change in effective diameter requires 
an axial deflection less than 5 mils per plate. Thus, for the total 
maximum axial compression of 0,15 inches (3.81 x 10 -3 meters), 
at least 30 plates would be required. This is more than can be 
fitted into the core length available for the symmetrically opposed 
configuration, which is less than 1.5 inches. It appears, in fact, 
that only about 12 plates can be used. Thus, the % AD can be no 
smaller than approximately 0. 26. e 
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TABLE JV 


Range of Parameter Variation for the Nonlinear 
Behftvioy of tho Symmetrically Opposed 
Toroidal Plates 


Axial Deflection, 6 


0 to 15 mils 

(0 to 3. 81 x 10 “ 4 m) 


Thickness, t 


3 to 7 mils 

(7.62 x 10 " 5 to 1.77 x 10“ 4 m) 


Pressure, p 


0 to 500 psi 

6 2 
(0 to 3. 445 x 10 newton/m ) 


Included Angle, ^ 60° to 180® 

Span, S (Approximate) 1/8 to 1/4 inch 

(3.17 x 10 " 3 to 6.35 x 10 " 3 m) 
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Effect of Plate Thickness. - The effect of variation in plate 
thickness on effective diameter for various included angles is shown 
in Figure 26. (Reference Run Nos. 3, 5, 6, 7, 8, 21, 24, 25, 26, 

27, 34 and 35 of Appendix D). 

It can be seen that plate thickness has practically no effect for 
the large included angles. The moderate effect shown for an angle of 
60® is somewhat academic for a bellows where constant effective 
diameter is the principle goal in view of the fact that its performance 
in this respect is not competitive with the larger included angles. 

Effect of Pressure Loading. - The effect of pressure loading 
itself on change in effective diameter for a 120® plate with a nominal 
axial deflection of 10 mils is shown in Figure 27. (Reference Run 
Nos. 21, 28 and 29 of Appendix D). It can be seen that, over the 
pressure range of 50 to 500 psi, the nonlinear analysis predicts 
very little effect on the change in effective diameter due to pressure 
alone. 

Effect of Included Angle. - The effect of included angle 
variation on the change in effective diameter is shown in Figure 28 
for two plate deflections. (Reference Run Nos. 2, 3, 15, 19, 21 
and 34 of Appendix D). This plot again shows that fact that a large 
included angle should be used for minimum change in effective diameter 
and that the variation is small between angles of 120® and 180®. 

Effect of Span. - The effect of span on effective diameter can 
be seen by comparing Run No. 2, having a span of 1/4 inch with 
Run No. 12, having a .1/8 inch span, in Table D-l of Appendix D. 
llqre, it can be seen that the change in effective diameter is 
approximately the same for equal axial deflection, all other factors 
constant. Thus, from a design point -of -view, if more plates of a 
smaller span can be fitted into the available length so that the 
deflection per plate is reduced, the expected percent change in 
effective diameter would be reduced according to the behavior indicated 
in Figure 25. Thus, it may be concluded that minimum span is 
desirable which, of course, is limited ultimately by manufacturing 
difficulties. 

'v Effect o f Boundary Conditions. - As discussed before, the 

solQl'NPfib^ased as a basis for the preceding discussion assumed edges 
which are both angularly and radially fixed. The wxigular fixity is 
provided automatically by symmetry for the opposed configuration. 
However, radial fixity depends upon the use of a sufficiently stiff 
ring welded between the bellows plates. The actual radial deflection 
of an opposed toroidal convolution plus practical stiffeners has been 
estimated as approximate y 1 mil. To determine the degree of effective 
diameter change, that might occur due to this actual radial deflection, 
a 1 mil radial direction was imposed at the boundary of both the 
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500 psi (3.445 x 10 newtons/m ) 

10 mils (2. 54 x 10 " 4 m) 

0.216 to 0,25 (.0055 to .0064 m) 
as noted on curve 

6 inch (. 1 52 m) 


0 = 60® 

S = 0. 26 in 
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S = 0. 216 in 
0 = 120 ° 

■pi= 180° 

S = 0. 25 in 




t ® 5 mils (1. 27 x 10 m) 

6 =10 mils (2. 54 x 10 ^ m) 

S = 0. 216 inches (. 0055 m) 

D =6 inches (. 152 m) 


1 1 i 1 
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Pressure, p, psi (newton/m x 1.45 x 10 ) 

Figure 27. Effect of Pressure Loading on Change in Effec 
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120° and 180° included angle plate (Run Nos. 13, 14, 30 and 31). 

The results of these runs are compared with the identical fixed cases 
(Runs No. 1, 3, 16 and 21) in Figure 29. It can be soen that the 
maximum variation of change in effective diameter occurs for the 
180° plate and is only 0.01 percent out of 0.2 percent. Thus, the 
error introduced by assuming a fixed boundary condition is negligible. 


Interpretation of Results. - The results of the parametric 
computer analysis summarized in the preceding paragraphs indicates 
that the smallest change in effective diameter that can be expected 
over the desired axial deflection and pressure range is approximately 
0.26 percent, within the overall design constraints imposed in Table II. 

It can be seen from Figures 25 and 27 that the nonlinearity 
or change in effective diameter is due entirely to the axial deflection 
of the plate and that its effects are fully present even at very low 
pressures. Two physical effects can combine to produce this type of 
nonlinearity: 

(a) A nonlinear force -deflection characteristic or spring 

rate at large deflection which would appear in the solution 
for combined loading as a change in effective diameter. 

(b) Distortion of the shape by the axial deflection or load 
so that the reaction forces produced by subsequent 
(or combined) pressure loading are distributed 
differently than with no axial deflection. 

A third possible effect, distortion of the plate shape under the pressure 
load itself can be seen by Figure 27 to produce little change in 
effective diameter. 

However, it can be seen from Run No, 9 of Table D-l. . 
(Appendix D) that the total axial reaction force due to a 10 millinch 
axial deflection alone is only 2.6 lbs for the 180° plate. Any non- 
linearity in spring force would be small compared with this total 
load. Run Nos. 1 and 3 for the same plate, on the other hand, show 
that the change in reaction force for this deflection at 500 psi pressure 
loading is 5. 8 lbs. Thus, it is obvious that this latter nonlinearity 
must be almost completely due to Effect (b) above rather than (a). 

The same conclusion can be drawn for the 120° plate from Run Nos. 

16, 21 and 23. 

Effect (b) could be produced simply by the plate tilt due to the 
axial deflection. This was investigated by Runs No. 16 and 17 of 
Table D-l where two identical cases were run with no axial deflection. 
They differed only in that the plate in Run No. 16 was tilted by 2.3 
degrees to simulate an axial load with, of course, no shape distortion. 
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of Effective Diameter, AD^ percent 


6 2 

p « 500 psi (3.445 x 10 newtons/m ) 



w.„ r . 29. Comparison of Chan.. In Effective nlameter B etween Change 
with a Fixed Edge and a Specified Edge Displacement.^ 
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Tho results for the two cases wore idontical and virtually linear 
over 500 psi. Thus, the distortion of interest is actual shape 
distortion away from the initial toroidal shape due to the bending 
stresses engendered by the axial deflection. 

Design Studies. - Several computer runs were devoted to an 
attempt to uncover design techniques for the opposed toroidal which 
would reduce the nonlinearity due to shape distortion under axial load. 

One general technique which was thought to be a possible improvement 
was to provide relatively compliant edges with a radially stiff span 
so that axial deflection is produced by plate tilt in order to minimize 
the plate shape distortion. 

One method for approaching the above condition is to vary the 
compliance of the plate across the span so that the edges are relatively 
much more compliant for bending deflections than the center portion of 
the span. This was done in Run Nos. 46 to 49 of Table D-l where 
the plate thickness and moduli were varied across the span in such a 
way that the circumferential stiffness was constant, but the meridional 
or radial bending stiffness was much greater at the center of the 
span that at the edges. For Run Nos. 46 and 47, the circumferential 
stiffness was equivalent to a 5 mils constant thickness plate but the 
bending stiffness near the center of the span was 9 times greater 
than that at the edges. For Run Nos. 48 and 49, the circumferential 
stiffness near the edge was equivalent to a 2. 5 mil plate and that near 
the center to a 5 mil plate but the central bending stiffness was 36 
times greater than that at the edges. The salient result was that 
the predicted change in effective area was exactly the same as that 
for a constant 5 mil bellows plate. The distribution of curvature change 
can be easily seen in Figure 30 where the angular deflection, (3 , 
is plotted across the span for these two cases and compared with a 
uniform thickness plate. As expected, the major effect was to 
concentrate the curvature change near the edges. However, it was 
so much greater in this region that any advantage from the stiffened 
central portion was offset. 

A final concept for a shape change. to approach the ideal dis- 
tortionless toroid was the use of a non -isotropic or radially corrugated 
plate. Here, corrugations in the radial direction instead of thickness 
variation or modulus variation as described above would be used to 
make the circumferential stiffness negligible compared with the bending 
stiffness near the compliant edges. This would have the effect of 
approaching a "rigid -body” rotation of the ideal toroidal shape and 
might reduce the effective area change. This type of configuration, 
however, cannot be treated using available nonlinear analytical techniques. 
Thus, it is recommended that a critical experiment be performed to 
check its effect on linearity. The basis for the design of a non- 
isotropic opposed toroidal plate is given in Appendix F. 
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Included Angle, <t> , Degrees 

Angular reflection Across the Span fo r the Ohlform Thickness a nd 
Compliant Edge Cases,, 





Another shape parameter that might be used to reduce the change 
in effective diameter is the selection of an optimum Initial shape so 
that the nonlinear effects duo to distortion compensate each other. 

This, in essence, means varying the initial curvature of the plate 
away from a perfect toroid in an optimum way. It might bo expected 
that variation of the shape so that the deflected shape would approach 
a perfect toroid might reduce the nonlinearity in reaction force under 
pressure load. This was tested by Run Nos, 42, 43 and 44, The 
results showed that the effective area change was identical to that of 
an initially undistorted plate. In fact, as might be expected, the 
nonlinearity is due to the shape change alone and minor changes of the 
starting point have no effect. 

A second curvature change aimed at producing a relatively stiff 
central portion similar to that produced by the variable thickness 
cases discussed in preceding paragraphs was tested by Run No. 45. 

Here a peaked plate with a relatively stiff central portion was produced 
by reducing the radius of curvature in the central region. As before, 
this concept produces no advantage in terms of linearity and has, as 
its main effect, the increasing of the peak stress in the plate. 


Nested Toroidal Plates 

The principal advantage of the nested toroidal bellows is that, 
although the individual plate behavior is not expected to improve over 
the opposed plate case, many more plates can be stacked into a bellows 
core having a given length. Thus, each plate is required to undergo 
a smaller plate deflection which has been shown in the preceding 
section to be the major variable affecting change in effective diameter. 

Another fundamental difference is that the convolution is not 
symmetric about a plane perpendicular to the axis as is the opposed 
case. Thus, the pressure load is expected to cause a rotation of the 
elements with a corresponding change in effective diameter even at 
zero axial deflection, 

A series of 50 computer runs were made to investigate the 
effects of varying a number of parameters on the performance of the 
nested toroid. Two basic configurations of importance, single and 
multiple convolution, were studied as defined in Figure 30 where fixe 
principal nomenclature is also defined. 

In order to minimize computer time, the single convolution 
(2 plates) case was used to determine the basic nonlinear behavior 
of the bellows. The multiple convolution case was then used for the 
more promising combinations of variables to study the effect of the 
end condition and select the most promising bellows plate and stiffener 
configuration. 
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It can bo seen from Figure 31 that, In general, a stiffener 
ring is provided at each plate joint and the effect of varying the 
stiffness of this plate was studied. In addition, all of the parameters 
studied for the opposed plate were also investigated for tho . 
plate. The only paramo^ fixed was the moan diameter of the bellows 


plate, 


D 53 
m 


6 inches (0,152 m). 


The linearity of the bellows was computed as expressed in 
TTrmaHnns (23) through (26). That is, the change in effective diameter 
^'proportional mthe change in axial reaction force. N p . between the 
linear and the nonlinear solution for a given loading. 

The results of the parametric study are summarized and discussed 
in the following paragraphs. The numerical results of the computer 
solutions are tabulated in Appendix E where sample stress distributions 
for the plates are also given. 

Effect of Axial Deflection and Included Angle. - The percentage 
change produced in the effective diameter between small pressure loads 
and a pressure loading of 500 psi (3. 445 x 10* newtons /m*) as a 
function of the axial plate deflection, 6 is shown in F \S* re32f ° r 6 
three included angles, (Reference Run Nos. 4, 5, 11, 12, 15 and 1 
of Appendix E). 

It can be seen that the change in effective diameter is again 
proportional to the axial deflection and has about the same overall 
sensitivity as the opposed toroid (related to individual plate deflection). 
The total change does not vary markedly with included angle m the 
range studied (also true of the opposed toroid xn this range). 
major effect of varying the included angle is to shift the curve dotard 
so that the deflection at which there is no change in effective diameter 
from the linear case occurs at a finite non -zero value. This ha 
some practical significance since, by selecting the optimum 0, the 
total range of change in absolute effective diameter over a range of 
deflections can be halved* 

Ettect of Plate Thicknes s. - The effect of variation in plate 
thickne ss on change in effective diameter can be seen in Figure .33 tor 
an included angle of 140 e . (Reference Run Nos. 17, 18, 21 and 22 
of Appendix E). It can be seen that, for this case, the increase in 
thickness from 5 to 10 mils results in a smaller absolute change in 
effective diameter over the pressure range. 

Effect of Span. - Although no specific cases were run to 
investigate the effect” of span alone, the behavior for one case can be 
seen by comparing Run Nos. 5 and 7 of Table E -2 where the span 
decreases from 1/4 inch to 1/8 inch, all other factors except plate 
thickness remaining constant. There was essentially no effect on 


68 








(b) Multiple Convolution 


- Span 

- Mean Radius 

- Plato Thicknoas 

- Toroid Radius 

- Included Angle 

- Stiffener Ring Length 

- Stiffener Ring 
Thickness 

- Axial Reaction Force 

- Axial Deflection 

- Radial Deflection 

- Number of 
Convolutions 

- Applied Pressure 


Figure 31. Schematic Diagrams and Definition of Nomenclature 
for the Nested Toroidal Bellows. 
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Figure 32. Effect of Axial Deflection on Cnange i n Effective 
Diameter for Single Convolutions with Various Included Angles. 




Change in Effective Diameter, AD , percent 






change in effective diameter which is the result expected from the 
single plate results of the opposed toroid study. This confirms the 
previous conclusion that, from a design point -of -view, minimum span 
consistent with feasibility of fabrication is desirable since more 
convolutions, in general, can be fitted into a given core length, 

Effoct of Edge Boundary Conditions. - The edge boundary 
conditions on the plate aro provided by the stiffness of the stiffener 
rings. Tho basic conditions used for reference in the single con- 
volution studies assumed radially and angularly fixed edge conditions. 
Here, the results are not markedly different from those of the opposed 
toroid study in terms of the range of change in effective diameter over 
a given range of single plate deflection. 

Several cases were devoted to checking the effect of relaxing 
the fixed boundary conditions. For example, the angularly fixed case 
for d> - 140 (Reference Runs 11 and 12 of Table E-2) is compared with 
the angularly free case (Reference Runs 19 and 20 of Table E-2) in 
Figure 34, both being radially fixed. This produced very little effect 
on the change in effective diameter* 

Also in Figure 34, the fixed case for <f> - 180° (Reference Runs 
4 5 of Table E-2) is compared with one having an imposed radial 

deflection of 1 mil (a reasonable value for typical stiffener rings) but 
remaining angularly fixed (Reference Runs 2 and 3 of Table E-2). “Hie 
total range of A D remains the saine for these cases, the main 

effect being a shi/t downward.^ This some 'practical value, as 

mentioned before, in that it reduces the^maximum absolute change m 
effective diameter experienced over the range of plate deflections. 

The main conclusion from these studies is that so long as 
reasonable radial stiffening is present, the stiffness of the edge 
boundary conditions does not have a major effect on the single 
convolution behavior. 

Ettect of Multiple Convolutions. - Because the nested toroid 
convolution is not symmetric, plate rotation must occur if the edge 
stiffeners are not sufficiently stiff to provide angular fixity. Thus, 
if the end fittings are very stiff so that the end boundary condition is 
essentially fixed some variation in plate behavior may be expected 
for succeeding convolutions away from the end fitting. The effect of 
multiple convolutions on the change in effective diameter was tested 
for a 160° toroidal plate by Run Nos. 23 through 29 of Table E-l. 

The results are plotted in Figure 35. Here, it may be seen that the 
two convolution case gives a result significantly different from the 
single convolution but that the difference becomes much smaller as 
succeeding convolutions are added. The major effect is to move the 
effective diameter versus convolution deflection curve down, thereby 
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Figure 35. Effect of Multiple Convolution on Change in 
Effective Diameter. 
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..etaining the same sensitivity of change in effective diameter with 
deflection but decreasing the range of absolute change. 

The specific dimensions for the plate and stiffener rings are 
shown in Figure E-l of Appendix E. The decay of the effect of the 
enu boundary conditions on the deflection of the plates can be seen 
graphically by comparing the plate stress distributions plotted in 
Figures E-2 through E-5. The stress distributions in the second and 
fourth plates are essentially identical and it may be expected that 
succeeding plates will show no differences until the end fitting on the 
other end of the core is approached. 

Thus, it may be concluded that 4 convolutions are sufficient 
for a good prediction of the effective diameter change in a multiple 
convolution bellows of this type where edge stiffener plates are typical 
of those of Figure E-l in terms of radial and angular stiffness. 

Interpretation of Results. - Comparing the basic single convolution 
behavior shown in Figure 32 with that for the opposed toroidal plate 
shown in Figure 25, it can be seen ti. at they are comparable from the 
standpoint of change in effective diameter over the desired pressure 
and deflection range. (It is important to note that the plate deflection 
scale of Figure 25 represents one -half the single convolution deflection 
for comparison with Figure 32). 

The use of a partial toroid (<£ - 120° or 140°) increases the 
sensitivity of change in effective diameter with deflection in each case. 
However, this increased sensitivity is accompanied by an adjustment 
in the zero change location for the nested plates. Thus, the net effect 
is to reduce the absolute change in effective diameter with deflection. 

The major difference between the opposed and the nested cases 
is that the latter are inherently capable of closer spacing, thus 
allowing more plates within a given core length. This reduces the 
required deflection per plate and, thus, the change in effective diameter. 
Since it appears that the 15 to 30 convolutions (30 to 60 plates) 
required to keep the change in effective diameter below 0. 1 percent 
for the various <b ' s can be physically fitted into the core length, this 
configuration, unlike the opposed toroid appears basically feasible 
for meeting the design goal. The analysis, of course, does not 
consider difficulties in fabrication which are much more critical for 
the nested configuration than for the opposed. These will be discussed 
in a later section. It also does not take into account the design trade- 
off associated with the variable number of plates for the various 
< p's that may be fitted into the bellows core. Several multiple - 
convolution solutions were devoted to this subject in order to select 
the most promising configuration as discussed in the following chapter. 


75 


SELECTION OF THE EXPERIMENTAL 
BELLOWS CONFIGURATIONS 

The objective of the analytical studios of the preceding sections 
of this report was the identification of those bellows configurations 
having the best promise for meeting the linearity goal, constancy of 
effective diameter within + 0. 1 percent across the specified deflection 
and pressure range. 

In fact, only the nested toroid appears theoretically capable 
of achieving this goal since sufficient convolutions can be used to 
bring the nonlinearity within this bound by reducing the individual 
plate deflection. 

However, the opposed toroidal configuration is theoretically 
capable of achieving constancy of effective diameter within the 0. 2 to 
0. 3 percent range with the number of convolutions possible. Since 
this is still much better than the present state-of-the-art for bellows 
seal applications, it was recommended as an alternative test bellows 
for this program. Also the opposed toroid is inherently easier to 
fabricate and was felt to have a considerably better probability of 
surviving this critical step of the development. Moreover, the 
radially corrugated version of tbe opposed toroid is felt to be promising 
for further reducing the effective diameter change, perhaps bringing 
it within the design goal. 

In s umma ry, therefore, the three most promising design con- 
figurations selected on the basis of the analytical studies for further 
experimental evaluation were: 

(a) Opposed Toroid 

(b) Radially Corrugated Opposed Toroid 

(c) Nested Toroid 

In addition to the nonlinear analysis discussed in the preceding, 
section, further studies and computer runs were devoted to the selection 
of the optimum configuration for each of these designs. These studies 
are discussed in the following paragraphs. 


Opposed Toroid 

The nonlinear analysis of single toroidal plate showed that the 
axial deflection per plate and the toroid included angle have the 
greatest effect on change in effective diameter. For the multi - 
convolution case, the number of convolutions that can be fitted into 
the desired core length becomes important. The maximum number of 
convolutions can be expressed approximately from the geometry of the 
bellows as, 
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where 


L sin (cft/2) 

n S (1 -cob 0/2) + i» 8 sin (0/2) 


L = core length 

L = stiffener length per convolution in excess of the 

8 toroidal plate height 

S = span 

0 = toroid included angle 

This is plotted in Figure 36 for a realistic span of 0. 25 inches and 
two stiffener lengths. In these plots, the fact that the number of 
convolutions must be integer is taken into account. Realizing that 
the individual plate deflection, 6, may be expressed by, 


A ~ 


A L 
2n 


the individual plate data of the preceding chapter can be modified to 
reflect the number of convolutions possible in the prescribed core 
length of 1.5 inches (0.038 m). This is plotted *n Figure 37 for a 
0. 25 inch (. 0065 m) span and a net stiffener ien ;th of 40 mils 
(1.02 x 10 "3 m ). The corresponding spring rates and combined stress 
are also shown in this figure. 

It can be seen that, although stress and spring rate are 
within the prescribed limits (150,000 psi and 800 lbs/in) across this 
range, the stress is minimized in the 120 to 160° range of included 
angle. The change in effective diameter increases monotonicaily in 
this range, making 120° the optimum included angle. This is a weak 
optimum, however, and the choice of included angle is affected from 
a design point -of -vie w by the fact that the 160° included angle allows 
a reduction in the number of convolutions from 8 to 6, making the 
bellows significantly easier and cheaper to fabricate. The increase 
in change of effective diameter is felt to be small enough so that 
the fabrication advantage is warranted. 

Thus, the 6 convolution 160 C toroidal plate is recommended 
as the best compromise for the opposed toroidal bellows. The 
change in effective diameter is nearly independent of plate thickness 
in the realistic 4 to 6 mil range. The recommended span is 0. 25 
inches and the mean radius is 3 inches. For the opposed toroid, 
the angular stiffness of the stiffener rings is not critical because of 
the symmetrical load placed on them. They must only be stiff 
enough to resist radial deflections. This is true with stiffeners having 
a net increase in convolution length of 40 mils as assumed for the 
data of Figure 37. 



Maximum Number of Convolutions within Prescribed Length 
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. For the corrugated opposed toroid, the 10 mil amplitude 
corrugation with a half wave-length of 50 mils described in Appendix F 
is recommended as providing significant stiffness ratios for experimental 

evaluation. 


Nested Toroid 

Selection of Stiffener Dimensions. - The results of the nested 
toroid analysis showed that, providing sufficient radial stiffening were 
present, the angular stiffness of the stiffening rings was not critical. 
Radial stiffness is provided by using the maximum cross-sectional area, 

From a design point -of -view, a maximum number of plates 
having mini mum span is desirable to minimize change in effective 
diameter. The minimum practical span was selected to be 0. 25 inches 
Since the maximum allowable span, including both toroidal plates and 
stiffeners, is 0.6 inches, the radial dimension or width, A R g of the 
stiffener rings should be approximately 0.15 inch. This was used as 
the nominal value for the recommended design. 

For the nested configuration, the stiffener ring length, L g , can 
be related to the total core length, L, by 


L (2n - 1) = L = 1.5 inches 

s 


(29) 


To ensure that adjacent plates still have an axial clearance, 6 c , after 
being compressively deflected through A! , the stiffener length must 
meet, 

l > — + 5 

s ~ cos (p 2 n c 


(30) 


where 

t = plate thickness 

<p = included angle 

Taking t * 6 * 5 mils, Equations (29) and (30) may be solved for 
various 0 to ^iv** L and n. These are tabulated m Table V, 

Selection of Included Angle. - Using the stiffener length tabulated 
in Table V and a width, AR , of l e " mils, the change in effective 
diameter was predicted for stacks of four convolutions, (Reference, 

Run Nos. 27 to 39 of Appendix E). The results of these computer 
solutions are plotted in Figure 38. 
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TABLE V 

Acceptable Combinations of Stiffener Length and 


Number of Convolutions for Various Included Angles 


Included Angle, 

<P 

degrees 


Stiffener Length, 
L s 

milli -inches 

Number of 

Convolutions, 

n 

25 

30 

30 

25 

35 

21 

40 

19 

55 

14 


milli -inches 


8 


10.7 


5 milli -inches 
5 milli -inches 
0. 15 inches 























It should be noted that a stack of four convolutions was used 
instead of the actual number corresponding to each included angle 
shown in Table V in order to save computer time and fit within the 
conveniently available memory capability. Also the preliminary 
multi -convolution solutions showed that four convolutions was sufficient 
to predict the basic behavior of the bellows. 

Also plotted in Figure 38 is a locus of the sing 1 ® 
displacement corresponding to the total deflection of 0 15 inches and 
the P selected maximum number of convolutions for each included angle. 

This convolution deflection is tabulated in Table Y* 

It can be readily seen by the points of intersection of the 
performance locus with the curves for various <p s that the goal 
0. 1 percent absolute change in A D is only met for plates having 
included angle in the range of 140 to 160 . 

The combined stress and spring rate for * h « B * decreases 

are clotted in Figure 39. Here it can be seen that the stress decreases 
^th^ncrea^ng ^ and that the spring rate, although increasing, is 
always well below the 800 lb/inch limitation. Thus, it appears tha 
the optimum choice for included angle is 160 * 

Selection of Plate Thickness. - In order to test t ^ e e ^‘ :t * ° £ a 
significant thickness increase on the performance of the bellows Rim 
Nos 40 through 42 of Appendix E were made using 8 mil plates. Th 
re suit s showed slight improvement of effective diameter change and 
stress but a nearly five -fold increase in spring rate. Since it is 
felt that this spri n g rate increase more than offsets any advan age 
gained 4 in con et&ncy of effective diameter and that the increased thmkn.ss 
Lay lead to increased fabrication difficulties, it was decided to use 
the 5 mil plate thickness. 

Fleets of a Welding Flange. - Later fabrication studies 
indicated the desirability of providing a lip or flange at ^ e ^ge of 
the toroidal plate for welding onto the stiffener ru *S 8, 
of this type of flange are shown, schematically m Figure . 

These were treated analytically by Run Nos. 43 through 50 of 
Appendix E. 

The circular flange of Figure 40(b) would not c °nve/ge to a 
solution having acceptable accuracy on the computer. Thus, a Uat 
flange mode^ls shoLi in Figure 40(c) was used and valid solutions 

were obtained* 

The general results indicated that there is no significant effect 
on change in effective diameter from this type of flange. Also, the 
spring rate is not affected. 
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Figure 38. Change in Effective Diameter for Various Included Angles 
with Acceptable Combinations of Stiffener Rings and Numbers of Convolutions. 
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Figure 40. Definition of Edge Welding Flanges Investigated Analytically, 


However, the peak stress is increased by a factor of over four 
(from 40,000 psi to 172,000 psi) putting it beyond the design limitation. 
Moreover, it is felt that this typo of configuration will also impose a 
stress concentration further increasing local stresses beyond the 
analytically computed values. 


Thus, it is recommended that the 
avoided if at all possible within the state 

use of this type of flange be 
-of -the -art in welding. 

In Summary. - It is recommended that the parameters used for 
the experimental nested toroidal bellows be as follows: 

included angle, 

0 

= 160° 

plate thickness, 

t 

= 5 mils 

number of convolutions, 

n 

= 19 

stiffener ring size, 

L . 

= 40 mils 


A®. 

= 150 mils 


86 


FABRICATION STUDIES 


A specific objective of the study was the experimental 
verification of the performance predicted by the analytical studies for 
the most promising bellows configurations. Thus, with the completion 
of the analytical work and selection of the opposed and nosted toroidal 
configurations, effort was turned toward the fabrication of special test 
bellows. 

The design of the bellows plate forming dies, welding fixtures, 
and chill rings (heat sinks), as well as the detailed design of t he 
bellows plates, their flanges, and the mating surfaces of the inter- 
connecting stiffener rings and end fittings was accomplished at the 
outset to produce optimum hardware implementing the basic designs 
resulting from the analytical work. 

Because each of the selected configurations described in the 
preceding sections requires the incorporation of stiffener rings into 
the bellows core to achieve the desired linearity, they presented a 
common fabrication problem, the joining of thin gage bellows plates 
with relatively heavy stiffener rings as shown in Figure 41. These 
configurations incorporate some modifications from the recommended 
configurations of the preceding section designed to facilitate fabrication. 
In each case, however, the modification was not felt to significantly 
affect the basic performance of the bellows. 

Several joining techniques were attempted to solve the problems 
associated with the plate to stiffener ring joint, but each was 
unsuccessful within the scope of effort possible in the program. 

This work is summarized in the following paragraphs in order to 
provide a basis for future development efforts. It is felt, at this 
time, that the stiffened toroidal bellows can be produced using welding 
techniques, but modified tooling in order to implement revised bellows 
plates and welding fixtures would be required to overcome the problems 
experienced. 


Inert Gas Welding 

The preferred method for joining the bellows elements was 
inei ; gis shielded -arc, non -consumable electrode fusion welding (TIG 
welcing). 

Assembly of bellows plates and stiffener rings in welding fixtures 
brought problems to light very early in the fabrication effort. A 
combination of fit -up requirements and optimum location of metal 
fusion resulted in unsatisfactory welds. Inner stiffener rings fabricated 
from annealed hot rolled Inconel 718 plate could not be produced flat 
enough to result in intimate contact for the full surface of the bellows 
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plate flanges within the stiffener ring centering recesses. Blow- 
through was experienced in those areas where the parts were not in 
contact. The use of stiffener rings machined from annealed forged 
rings in order to effect radial grain flow and reduce the effects of 
the orientated grain structure of tho rolled material, was indicated. 

Bellows plate weld flanges were designed with what was con- 
sidered the maximum practical radial dimension since the desirable 
characteristics of the toroidal type convolutions were predicated on 
close control of the geometry at the juncture of the plates and stiffener 
rings. The geometry of the stiffener ring as shown in the detail 
of Figure 41 was chosen so as to produce optimum wetting and fusion 
of the two elements in order to produce a weld as shown in Figure 42, 

With 4 mil (1 x 10 ” 4 m) toroid plate thickness and 10 mil 
(2.5 x 10" 4 m) radial flange width, fixturing which would adequately 
centralize and bottom the plate within its recess and still allow 
adequate clearance between the welding electrode and the fixture 
elements, required great attention to detail. 

Numerous fixture changes and rebuilds were made and the 
effects of increased ring recess depth were evaluated by welding 
elements incorporating these various changes. All were considered ineffective 
as solutions to the problem. A continuous joint of the quality depicted 
in the photomicrograph cf Figure 42 could not be consistently produced 
throughout 360°. 

In the light of the lack of success as recorded above, it was 
necessary to consider alternate means of fabrication by TIG w siding 
techniques and the possibility of alternate welding systems. 

A design change compatible with the TIG system wae implemented 
for trial. This involved the addition of an integral lip on the stiffener 
ring, so located that it might be mechanically rolled over the bellows 
plate, flange, retaining the flange in its stiffener ring recess. Figures 
43(a), (b) and (c) illustrate the original and the changed design. 

This innovation gave great promise on the first trial weld. 

One plate was satisfactorily welded throughout the 360° pass. This 
plate was an opposed type toroidal segment plate and was welded to 
its . 200 inch (5 x 10" 3 m) thick inner stiffener. The assembly sequence 
required that another similar plate be welded to the opposite side of 
the same stiffener ring. In the process of rolling down the retainer 
lip on this plate, the clamping action or lip -to -flange overlap was 
minimal over part of the circumference of the flange -stiffener ring 
interface. This was caused by out -of -roundness developed in the 
stiffener ring by localized heating during welding of the first plate. 

The actual clamping was not sufficient to retain the plate in position 
during this second welding operation. The same end result developed 
on parts having both plates lip -clamped in place prior to welding of 










cither plate. Changer* in flxturing wore no) effective in the reduction 
of this problem. 


Sevoral trial wolds wore produced with variations in iixturing 
and a tiff oner ring design. Perfect 360* wolds were not obtained 1 Mth 
nny degree of consistency in the production of accoptoblo innor stiffen! r 
ring and plate assemblies, 


Electron Beam Welding 

As an alternate welding system, Electron Beam Welding was 
selected as most promising for trial. Segmented toroidal bellows 
plates and inner stiffener rings for both the nested toroidal assembly 
and opposed plate assembly were made available to a subcontractor 
for production of specimen welds. Temporary fixturxng was fabricated 
bv the EB Welding contractor and trial welds were produced as 
depicted typically in Figure 44. As in TIG welding, parts . 
inadequate with fusion not occurring at the radius joining the bellows 
plate flange and the arcurate form of the semi -toroid. However, the 
experiment indicated an excellent chance of success with proper tooling 
and improvement in the form of the stiffener -ring recess. The sub- 
contractor 1 s estimate of the funding required to fabricate proper 
permanent fixturing, setting up the EB equipment and producing the 
actual welds required for bellows assemblies was beyond the funding 
available in the program at this point in time. 


Brazing 

Another joining system, brazing, which had been considered 
earlier in the program and rejected because of anticipated lower joint 
strength as compared to the welding proce-sses was at this point in 
the program given serious reconsideration as it was agreed that 
evaluation of bellows stability performance at less than the original 
test parameters would provide a contribution to the program. 

General brazing experiments were performed by a source 
specializing in High Temperature Vacuum Brazing. The brazing alloy 
recommended as most suitable for production of a joint having maximum 
strength and compatability with Inconel 718 alloy and its final heat 
treating temperatures was 82-18 gold -nickel. The first process 
specimens consisted of an opposed type toroid brazed to an inner 
stiffner ring. Again, the trial of an alternate joining system produced 
promising results. Figure 45 is a photomicrograph of an early specimen. 
An excess of brazing alloy occurred at the critical point cf the joint, 
the juncture of the outer diameter of the inner stiffener ring and the 
arcurate form of the semi -toroid. Several joints were produced xn 

which the control of brazing alloy in the critical area was given special 
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attention, Figure 46 shown the progress made In thin direction, 

Figured 46 and 46 dhow crass -sections of the sub-aBScmbly of 
an Inner stiffener ring and two semi -toroid plates. The attachment 
of the required number of such sub -assemblies to intermediate outer 
stiffener rings as shown in Figure 47 was the next step In the 
fabrication of teat bellows, The dotailod parts wore delivered to the 
brazing source for production of toroid -inner stiffener ring sub- 
assemblies which it was expected would be TIG welded to outer 
stiffener rings, TIG welds in the outer joints were expectod to present 
no problem since they were designed geometrically similar to conventional 
Sealol bellows end fitting welds. Reference to the outer joint in 
Figure 43(c) will further explain the geometry of these welds where 
the weld bead joining the parts are formed by the fusion of elements 
having a reasonable ratio of cross-section. This ratio is in the range 
of two or three to one. 

Sub -assemblies were brazed and tested for leak -tightness. 

These units required the repair of minor leaks in the brazed joint. 

This reheat, while it was performed in a controlled atmosphere 
furnace, produced an oxide layer on the surfaces of the parts. This 
oxidized material defied removal by bake -out in the vendor's vacuum 
furnace, nor could it be removed by a reverse plating process which 
produced unacceptable reduction in toroid plate thickness. 

Although the original intent had been to TIG weld at the plate 
O. D. ' s the brazing vendor suggested brazing all joints in each 
assembly simultaneously as a means of avoiding furnacing the sub- 
assembly units more than once. Trial of this suggestion would have 
required the production of additional detail parts and the purchase of 
considerable labor for pre -brazing preparation and expensive furnace 
time which could not be funded by the program at this point in time. 


Surface Contamination 

In addition to the problems associated with positioning of the 
parts and accurate fusion location control as described in the preceding 
sections, the quality of the fused metal was less than fully acceptable. 
This proved to be due to contamination traceable to tenacious surface 
oxides present on the 718 Inconel 4 mil plates fabricated from strip 
as received from the reroller supplier. Cleaning, pickling and 
storage in evacuated containers would not provide the program with 
contamination free plates. Experimentation did develop a successful 
electro -chemical reverse plating means for this purpose suitable for 
use on Bmall numbers of plates. This process also resulted in a 
means of plate thickness contro* through the removal of surface 
material. 
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BRAZED 


Figure 46 Braze Joint Specimen With Acceptiblc 
Amount Of Bra/e Alloy 
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BRAZE JOINT 



OPERATION 

1. Weld Two I. D. Stiffeners 

2. Braze I. D. Stiffeners ( Halfs ) 

3. Weld or Braze Plate to O. D. Stiffener 

NOT TO SCALE 


Figure 48, Improved Weld Construction 


DESIGN AND FABRICATION OF BELLOWS 
TEST APPARATUS 


A toot apparatus was designed and constructed for tho 
of verifying theanalytical results under cryogenic conditions. Although 
apparatus was not used for bellows testing in this program, it 
is described briefly in this section for reference purposes. 


Specifications for the Test Apparatus 

The basic purpose of the apparatus is, of course, the measure- 
ment of the change in effective diameter for a test bellows over the 
derating ra ng e o § f interest. Since the effective diameter or areas and 
pressures of interest combine to produce large pressure forces m 
comparison with the change due to nonlinearity, the approach adopted 
was to balance the pressure force on the bellows against that on a 
fixed, known area built into the apparatus and measure the difference. 
Since the expected variation in effective diameter is quite small 
(0 1 to 0 3 percent), the apparatus must be quite sensitive to variati 
in pressure force and must introduce very little friction force at its 

seals. 

The pressure range applied to the bellows was specified to be 
0 to 500 psi (345 x 10® newtons /m 2 ) applied internally and externally. 
The test bellows are to be deflected axially over a range of 0 to 
0.150 inches (3.8 x 10 " 3 m). The fluid acting on the bellows will 
be both liquid and gaseous nitrogen. 

In addition to the static tests above, a fatipe test consisting 
of cycling the test pressure with amplitudes of 0 to 500 P® 1 1*45 x 1 
newtons/m 2 ) with a constant applied bellows deflection of 0.150 inc e 
(3.8 x 10 " 3 m) was specified. 


Description of the Apparatus 

The apparatus, which was built and proof tested in parallel with 
the development of fabrication techniques is shown in ^e photegraphs 
of Figure 49. Figure 50 shows an overall schematic diagram of the 
test fixture. The description in the following paragraphs is keyed 
to the element numbers shown in this figure. 

The overall structure of Figure 50 consists of the main support 
frame, base plate, and rods supporting a yoke which provides means 
for adjusting the axial position of a load cell relative to a basic 
pressure vessel in order to obtain load readings at various bellows 
test lengths and pressures. 
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The prcoouro vobdc) is insulated from the frame base by 
Transits (1) which supports a steel pressure structure ar base plate (2), 
Directly attached to this structure are a ptatic separable outer pressure 
cylinder (3) and a multi -olemont static inner pressure cylinder (4) and 
the lower element of a Bellows Test Specimen (5). Fill ports (6) 
for the introduction of the cryogenic LN 2 to the inner and outer surfaces 
of the test bellows are provided in the base plate. 

Flow spreaders (7) or baffles will direct the flow of the LN 2 
to minimize the rate of localized cooling in the test bellows chamber. 

Fluid level sensors (8) were designed into the outer cylinder 
wall to indicate the level of cryogenic liquid in the chamber outward 
of the Test Specimen. The sensors are upper and lower elements 
of an ultrasonic system indicating the range of the LN 2 liquid level 
and are located to indicate complete immersion of the test specimen 
in the LN 2 and provide assurance that the upper L 1 N 2 level does not 
exceed its design limit. Control of the LN 2 level in the chamber inward 
of the test specimen is accomplished by a similar set of sensors (9) 
built into the base (2). In each case, thermal isolation of the external 
(16) and internal (17) diaphragms is maximized by insulating elements 
(10) and (11). Insulation (12) and (13) is also designed to assist in 
the isolation of the upper section of the apparatus. These elements 
are attached to cylinders (14) and (15) which transmit the forces 
developed by the test pressures on the specimen effective areas to 
the strain -gage load cell (25) which is adjustable in vertical position 
to effect deflection of the Bellows Test Specimen through elements 
21, 20, 19, 18, 15 and 14. All elements of the pressure vessel are 
thus to be cooled and temperatures of the apparatus stabilized. 

Covering (24) was designed to provide a major means of reduction 
of heat input to the cold bellows test chamber. 

Figure 51 shows an enlarged schematic diagram of the 
differential force sensing elements of the apparatus. The fixed 
reference areas for measurement of differential forces are defined by 
sealed cylindrical surfaces that allow at least 0.150 inch (3.8 x 10 “ 3 m) 
relative axial movement between surfaces without friction loss. Self- 
compensating rolling diaphragms (16) and (17) were selected as seals 
having the desired characteristics. These elements, to operate properly, 
must be isolated thermally from the cryogenic temperature in the 
chamber housing. This is accomplished by means of the insulating 
elements described in the preceding paragraphs. 


Operation 

In operation the upper section of the pressurized side of the 
vessel relative to the Test Bellow b will be filled with Gaseous Helium 
at test pressure to minimize loss of LN 2 and for further protection of 
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INNER HOUSING 
( FIGID ) 


) 

Figure 51. 


OUTER HOUSING 
( RIGID ) 


Differential Force Sensor 





the rolling diaphragm seals. Gaseous nitrogen will be vented to the 
atmosphere from the low pressure side of the vessel, Thcso conditions 
will, of course, be reversed as the test pressure conditions are 
reversed relative to the bellows. 

Test Pressure measurements are intended to be read on a 
precision pressure gage sensing the pressure of the GHe above the LNg 
on the high pressure side of the test bellows. Loads are to be read 
on a digital indicator of the output signal from the strain gage load 
cell, Refering to Figure 50* Test Bellows deflection is determined 
relative to a Reference Height Pin (23) representing the compressed 
height of the bellows. Adjustment of the bellows height on test will 
be made by adjustment of the position of the load cell and actual 
measurements of the adjustments of the position sensing element 
(19) relative to the Reference Height Pin (23). 

The apparatus was designed to perform Bellows Fatigue Tests 
by cycling the pressure applied to the LN£ external to the bellows 
over a range of 0 to 500 psi (345 x 10° newtons/m 2 ) at a rate of 2 
cycles per minute and with the Test Specimen in its compressed 
position. The GH3 blanket above the LNz fluid will be the pressure 
control medium handled through suitable valving and timing devices. 


Testing 

The Test equipment was constructed during the period of bellows 
fabrication and proof tested. The Apparatus was subjected to 750 psi 
(518 x 10® newton/m 2 ) water pressure at room temperature applied to 
the inner and outer pressure sections with a solid cylindrical element 
installed in place of a Test Bellows to separate the chambers. The 
only problem encountered concerned the rolling diaphragms (16) and 
(17). Clamping of the flanges of these parts proved inadequate as 
indicated in Rgure 52 which is a photograph of a diaphragm which 
pulled loose from the clamping surfaces. Roughening of all clamping 
surfaces for each of the diaphragms shown in the partial view of 
Figure 51 for the upper section of the apparatus proved to be the 
solution to this problem. After this measure, the apparatus 
successfully passed the Proof Pressure Test. 

Calibration Tests were not performed on this equipment as 
fabrication of test bellows was not successful. Verification of diaphragm 
areas was to be made against a telescoping cylinder having a sealed 
sliding diameter comparable to the designed areas of the diaphragms. 
Cylinder seal frictional hysteresis was to be obtained as the difference 
in forces in compression and extension and removed to correct the 
net forces, to be used in calibration of this element. Tests performed 
at room and LN£ temperatures would indicate temperature sensitivity 
of the apparatus and need for the addition of temperature sensing. 
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Figure 52. Subsequent to Pulling Loose 
From the damning Surfaces 
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CONCLUSIONS 


The general overall conclusion that may be drawn from the 
analytical studies performed in this program is that for the range of 
bellows pressuro, deflection, and size specified for this application, 
only toroidal bellows eloments appear promising for achieving 
significant improvements in holding constant effective diamotor. 

Opposed toroidal plates appear capable of holding their effective 
diameter constant within 0.2 to 0.3 percent. Nested toroidal plates, 
because of their capability of providing more convolutions within the 
envelope appear capable of meeting the 0, 1 percent goal. Even the 
opposed toroidal performance, should it be attained, would be 
significantly better than the current state of the art in seal bellows. 
Both configurations require a core made up of alternating flexible 
bellows plates and stiffener rings. 

The fabrication studies pointed out the difficulty of fabricating 
this type of structure in the desired size range. No successful 
bellows was produced within the scope of work allowed by the contract, 
although several individual welded joints were successfully produced. 

The nested toroid, as was expected, proved much more difficult to 
fabricate. The main conclusion reached as a result of the fabrication 
studies is that the stiffened opposed toroid can be fabricated with 
further development of the configuration near the bellows plate -stiffener 
ring interface and improved fixturing. Considerable attention must be 
given to improvement of the surface quality of the Inconel 718 foil 
material and the use of light oxide removal procedures in order to 
achieve good welds. 

It has been concluded as a result of the test development work 
that the pressure forces on the bellows in the desired range can be 
measured accurately enough to establish performance capability under 
cryogenic conditions. An apparatus for this purpose has been designed 
and built. 

The systematic evaluation of the range of possible bellows plate 
configurations for this application resulted in several detailed con- 
clusions which are summarized below for completeness: 

(a) The performance predictions for the conventional 
corrugated plate bellows showed that it was 
extremely marginal in this application for meeting 
the nominal performance goals (spring rate and pressure 
resistance). Thus, it was not investigated in detail 
for non linear effects, and it was felt that corrugated 
plate bellows do not offer promise for achieving the 
small effective diameter change goal. 


The flat plate bellows suffers from the same factors 
that make the corrugated plate unsuitable in the 
specified size range, 

A wide variety of S -shaped plate configurations wore 
analyzed and shown to have nominal performance 
characteristics comparable to, but generally inferior 
to, corrugated and flat plate bellows. 

Analysis of orthotropic plates (corrugations in a 
radial direction) showed that this feature has some 
promise, particularly on toroidal plates, for improving 
their linearity. 


recommendation 


The bellows designs developed In this program are expected 
to significantly improv the state of tho art in achieving small change 
in offoctivo diameter fo* cryogenic seals. Accordingly* it is 
recommended that this work be used as a basis for continuing devolop 
mont toward this end. 

Specifically, the next phase of the work should be aimed at 
developing improvements in the techniques used to weld flexible 
plates and stiffener rings having dissimilar sectional thicknesses. 

This will require modifications in both bellows plate and stiffener 
ring configuration near the interface. Also, improved weld fixturing 
will be required. 
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APPENDIX A 


DESCRIPTION OF COMPUTER PROGRAM "SEALOL" 


Computer program SEALOL prodiets the stresses and deflections 
of a pair of metal bellows plates under combined axial force loading 
and prossuro differential. It uses the linear shell equations and, thus, 
is limited to small deflections if an accurate solution is required. i*or 
the present study, SEALOL was used to conduct the screening analysis 
of corrugated and partial toroid plate shapes in order to establish 
their basic performance capabilities. 

In this appendix, the capabilities, the programming basis, and 
a listing of SEALOL are provided. 

Since NONLIN (described in Appendix B) is also capable of 
providing a linear range solution as well as a nonlinear one, a brief 
user -oriented comparison is given between the two programs in 
Appendix C. 


Capabilities 

SEALOL can treat bellows convolutions consisting of the various 
plate shapeB shown in Figure A-l as well as combinations of these 
cases. In the case of toroidal arc segments, up to four different 
segments having different radii can be treated. 

The basic equations allow only constant plate thickness across 
the span. Stiffening rings can be represented only by imposing simple 
boundary conditions on the plate, (e. g. radially or rotationally fixed). 

Both differential pressure loading and axial deflection can be 
treated. 


Governing Equations 

The program uses the classical firBt order shell equations 
reduced to two linear second order differential equations in terms of 
angular rotation, p , and horiaontal force, rH. This formulation is 
similar to the ReiBsner -Meissner transformation. 

The governing equations can be written, 
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(a) Conical Plates 



(b) Pairs of Conical Plates 



(c) Flat Plates 



(d) Toroidal .Arc Segments (inside weld) 
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The nomenclature for these equations is defined in Table A-l 
and in Figure A -2. 

The various possible prescribed boundary conditions for the 
plate edges are: 

(a) Prescribed Rotation or Moment 
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(b) Prescribed Radial Displacement or Force. 
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table a-i 


Nomenclature 


A 

P 

6 

D 

E 

h 

H 

P 

V 


E3 

3 

Q 

3 

a 

a 

2 


S 


</> 


r i 

r 

o 

r, z 




u 

V 


( )' 

( )" 


1 

Eh 

rotation of middle plane 
axial dieplaeomont 
Eh 3 /12 (1-P 2 ) 

Young* 8 modulus 
thickness of plate 
radial force 
pressure load 
Poisson’s ratio 
rH 

inner radius 
outer radius 
cylindrical coordinates 
axial load 

radial displacement 
axial force 

_JL ( ) 

0r ' ' 

n2 . . 



tfj a r o H Q or r^H^ at r 3 r Q or r^ 


Method of Solution 

The general computational method consists of the following steps: 

(a) The computer program uses finite difference equivalents 
for the governing differential equations (central difference 
in the interior and forward or backward differences at 
the edges). 

(b) The program determines a set of influence coefficients 
for each plate to represent the various deflections as 
functions of the imposed loadings and edge reaction 
forces, 

(c) These coefficients are then used to formulate a set of 
simultaneous linear equations for each element of the 
bellows, 

(d) This set of equations is solved for the unknown 
reaction forces taking into account the prescribed 
deflections or loads. 

(e) The resulting known edge conditions can then be used 
to compute the displacements and stresses at any point 
in each bellows plate. 

The computer program consists of the main program plus six 
subroutines. A logic flow diagram, defining the functions of the 
subroutines is provided in Figure A -3. The FORTRAN II listing of the 
program is given in the following section. 
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Road Input Data 


Transfer Design Input to 
Computer Input 


Compute Shape Function 
Using Subroutine ZFR 


Compute Geometric Parameter 
Using Subroutine AFR 

r z. 

Set Up Matrix Equation With 
Subroutine YBMAT and Solve 
Using Subroutine SOLV 

J — 

Compute Derivatives Needed to 
Get Stresses Using Subroutine PER 


t 



Compute Displacement Using 
Subroutine INTEG 




l - 


Compute Stresses and Influence Coefficients 


r , 

Solve for Final Solution and Print Out 


Fipure A -3, Logic Flow Diagram of Program SEALOL. 
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NWag 

NGe«a 
NGE*3 
NWCl* NHC?Bl 
NHC1 ♦ NBC?b2 


MAIN PROBl.FM 

NR® 1 PROBLEM NORMALLY CONTINUF 

CALCULATION for FINAL RESULT OF PLATE A only 
ELFmEnT a UF SEALoL BELLOW PLATE 
ELEMENT P OF SEAOOL RtLLOW PLATE 
FLAT HELLUw FOR TESTING 

wiVfeu MOMtNl ♦ OlVEN HORIZONTAL FOPOt 
wiVEN ANGULAR OISP.* GIVEN RADIAL DISP. 

NHCl* MHC?®3 DiVfcN ANGULAR UlsP** GIVEN HORIZ* FORCE 
NH" 1 , imBC 2 bA uIVEN MOMENT* OlyEN RADIAL Dls p . 

BClt HC3“ HOME "i I OR ANOULAR UIS p » 

BC2* RCA® HORlZ. FORCE OR RaOIaL OISP. 

NBCl, MCI, HC2*» INSIDE BOUNDARY 

NHC2* BC3* HCA® OUTSIDE BOUNDARY 

NCOnTbi* only for given boundary condition 

NCONtbZi CALCULATE the INFLUENCE COEFFICIENTS 

NCONTb.I* calculate the complete SOLUTIONS FOR TWO PLAtE 

STOP 111 ;NRC tKROR 

STOP AAA SINGULAR MATRIX 

SP IS SPAN 

SENSE SWITCH 1 U»PRINT OUT N£W Y MATRIX, U®CONTINUg 
SENSE SWITCH 3 UbcIRCULAH aHC, UbSEALOL BELLOW 
DELT® CORRECTION DUE To THIcNnFSS 
CPTsNOPMALIZEu pitch 
load»i Inside support 

LOAD®? OUTSiO.E SUPPORT 
LOAD®3 COMPLEILY CLOSED BELLOW 

DIMENSION C(lu*10)* 2(150) , DFL ( 150) » U(150). SRO<150>* SCDU501* 
1 SRBUSO) ,SCB(l50) »Ca (10,10) »CB (10,10) ♦ 

1 R(iso), DZ(iPO), UPZ ( ISO)* AU150) » A2 . „ 

j AA(150), AS ( 150 ) , Aft ( 150 ) * A7<150), F(150>* G(150). 

I r>F ( 15u ) * DG ( 150 ) * ZNEW(IOO)* 

j pMOMT(lOO) , CMOMT(IOO), FH(lOO), COE1(10*10>* CUEZ(10*10> 

ltCOE(lOtlO) , PR U0,2) ,E«(1V),B0(A) , 
lOSP (25,2*2) , US (50,2,2) ,H0t (SO, 2,2) , 

|RDs(5n,2,P) , CDs (50,2,2) *««S (*0*2*21 *CBs(50,2*2) » 

IRMMT (50*2,2) »tMHT (50,2,2) »FHS(50,2,2) 

COMMON R, DZ, DU2, F, 0, OF* 0G* CVP* Al* A2* A3* AA- AS* Aft. A7» 
1 XNIJ* U* A* CONS I 

3 READ INPUT Tape 5*13* ( (LA(I*D) *I®1*7> *J®1*A) 

READ INPUT TAPE 5*13* t (CB(I*J) * I® 1 *T) *J®1*A) 

®EAO INPUT TAPE 5* 1 0 * SPI * W ATA *RATB 
READ INPUT TAPt 5*10 *ES*XnU*SV 
A READ INPUT TAPt 5*10 *SP*RI*TH,dELT*CPT 

READ INPUT TARE 5*)2*N,NBCl,NBC2*NC0NT*L0AD 
00 TO (5*BOO*BUO) *MCONT 
5 READ INPUT TAPE 5* 10»5CI*5C2*BC3*BC4 
SoO READ INPUT tare b.BOl.OL.PRLD,LCONT.lRUN 
10 FORMAT (6E12*5) 

8ol FORMAT ( 2 E 12,5*215) 

12 FORMAT (1015) 

1 3 format ( 7Fio, 5) 

PLS®-1,00 

PIL®0.00 

NGE®1 

L«»l 
01 0 NR®1 

WRITE OUTPUT tape 6,812* IHun 
8)2 FORMAT (///16H A**RUN NUMHER®IS //> 
ft T®Th 
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PL»PLS 

Pl*PlL 

E«f8 

8i) TO 

1& DO 16 J*l«3 
DO 16 1*1*7 
)fe C ( I * J> "CA ( 1 * J> *PP 
DO S16 J*4*b 
DO 516 1*1*7 
Bl6 C(X.J)*CAn*J> 

PTA*(SPT"CPT>*BP«*RATA 

DO TO SO 

18 DO 19 .1*1*3 
DO 19 1*1.7 

19 CU»J)*C6(!*J:*^P 
DO 519 J»4*5 

DO 519 1*1*7 
619 C( J*J)*CB.(1« •' 

PTB*- (5PT-CP7 ) # 6P*RATB 
WHITE OUTPUT 6*21 

?0 PT»CPT*SP 
RO*RI*SP 

A*1 ./ (E*T) , 

01*12. *U.-XNU**2.) 

DbE*(T**3*)/01 
D IS STIFFNESS OF PLAT| 

WRITE OUTPUT IMP* 6* 923 
80 TO (901 *90«! *903) ♦ LOAD 
9ni write output **924*pl*pI 

P*J»L ' 

Go TO <*2Q 

9pS WRITE OUTPUT lM»-t 6*925*PL*Pi 
P«PL»PI* (R0**S"K1**2) *3* 1*159 
80 TO 920 

9o3 WHITE OUTPUT 1*PE 6.926*PL.Pi 
P»P|L-PI* (RO*.*S ) *3 • 14159 
60 TO 920 


UOAUINO CONDITIONS) 

SUPPORTED INSIDE AND AXIAL LOAD 

INTERNAL PRESSURE LOA0*E15.5> 
SUPPORTED OUTSIDE AND AXIAL LOAD 

jNTFRNaL PRESSURE L0AD*E15,&) 
CUOSf.D BELLO* AMD AXIAL LOAD “ E18 * S ' 

INTERNAL PRESSURE -E15.5) 


■E15.5/ 


*?3 FORMAT i//aOH _^*!!5 K c VK^^2n AXIAL load -eis.b/ 

)?4 FORMAT (/ASH 
146H 

>?5 FORMAT (/47H 
14TH 

i?t FORMAT (/39H 
139H 
9?l) PL«R 

WRITE OUTPUT TAPE 6*22 

WRITE OUTPUT TMPfc 6*23*R0*RI*SP*PT*TH 

WRITE OUTPUT TMPE 6.24»PL*E* XN(J*D*A 

J*N«*2 

„ "»<•«" INC.. 

» fSSSSli|S M %o:FB;sU» H “ellfMeiSS 1 FlTCH.re.M13H 
H*F«J!??S5’“t5;!lEi....*» «»».*•»" «NN«.».e.»H ceie.s.s 

le’ropeM'zTH 51 nuhber cr division, n.s.ibi 

50 N6ES*N0E , 

60 TO (52*53»TU*53) *N8t 

52 Pt*PtA*DELt 
80 TO 60 

53 PTapTB-OELT 

nG£«2 . _ . 

60 CALL SEAL0L(C.Hl*H0.PT.N6t) 

NQE-N6ES 

CALL 7FR(C*N*h*P*/) 
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BQ TO HO 

XN«N-? 

HbSP/XO 


• WHITE OUTPUT IhHE 6*73 
73 FORMAT ( 32H0 ^ U a I MELLOW 

DO 76 l ■ l » J 
XIb I 

m I )■«!♦( XI"H« > "H 

ZU)oO.O 
D2< I)»0»0 
76 DO7<n»0.n 

aO CALL AF« C R • 2 AMU * 0 * M 


PLATF FOH TESTING) 


) 


Hi KONT-t 

00 TO (»A* 9b >♦•*)< 
A* 00 TO <99.H5*8b> * 
Hb BClaO. 

WC2*0. 

riC3»0. 


NCOUT 


BCAsO. 

P»0. 

0o"tO<H 6*P7*0«*89*9|)*91) *KONT 
Qb BC I b 1 * 

00 TO 99 
fi7 BC2 b1. 

00 TO 99 
P8 BC3»I. 

00 TO 99 
ft9 BCAbI , 

00 TO 99 

90 P»l. 

00 TO 99 

9 1 PI»1. 

00 TO 99 
9 b BCla-HCl 
HC?b-bC2 
BC3s-Btf3 
BCAa-ACA 


PBPL 

PIopIL 

99 WHITE OUTPUT TAPE 6*25 
WHITE OUTPUT TMHfe 6*28 
00 TO (30*31*3**33) »NHCi 

30 WHITE OUTPUT TAPE b*3b*bCl*BCg 

31 WHITE OUTPUT TAPE 6«36*bCltf»C2 

00 TO AO _ 

32 WHITE OUTPUT T«t*E 6*37 *oCl *002 

00 TO 40 , 

33 WHITE OUTPUT T*PE 6*3H*OCi*BC2 


GO TO AO 
pb FORMAT (22H 
90 FORMAT (8IH 
35 Format ( 9H 
3t> FORMAT (16M 
3? FORMAT I lbH 
30 FORMAT (9M 


oUUNUAPY CUNUIT10NS) 

a( InSIUE BOUNDARY) _ .... 

Mot'ti ,, TBEi3»5*i9M horizontal forc£bEi3 
angular dism»*ei3.5*19h radial D!Sp**ei3*»| 
Angular DJSM.*E13.B*I9H HORIZONTAL F0RCEbL 13*6) 
MO'tN t«t i 3,6* lbH RADIAL DIs^» b EI3,6) 

AO ’white’ output Tape 0**f 

GO TO (330*331 *332*333) *NbC2 
A2 FORMAT ( 19H0 OUTSIDE BOUNDARY) 

30 WRITE OUTPUT UPi 6*3S»b63*bCA 
60 TO 100 

31 WHITE OUTPUT TAPE fc*36*bC3»BCA 
00 TO 100 

132 WHITE OUTPUT TAPE 6*37*BC3*BCa 
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313 

lpy 

6?9 


l?o 


00 TO TOO , 

WHITE OUTPUT l*Rt h»3h«or..4.MU4 

on to iou 
sa VF» sy 

WHITE OUTPUT TaRR Atb39*P*PI , „ _ 

FORMAT </l4H MAIai, UnAU»fel3.B,2lH INTERNAL PRESSURE »U3. 5' 
CaLL VRMAT (f tC*»N*h«KNUt a tUtCtPCl *NhC?»SV *P«PT »RCt ♦MC?*BL3 iHCA i 
SVaSAVE 

CALCULATION OR VERTICAL OfcFI fcCTlON 
Call HT 6 flfF*Mvi«»UEL> 

JaN/2 

T«RL A, 122 
UR£ Atl23 

TAKE a» 12 A* < w ( 2 *! ) *DFL( I > * I * 2 * J) 


HO 


112 

133 


Take At 133 

| ARE Atl 2 A» <H( 2 *I) tZNEW(I) tl« 
/U) AfTtR OtFO«M A TION> 

H ZNEW 

ZNtW R 

VERTICAL IlfcFLEC • ION) 


nEL 


znEw) 


dpl> 


ZNEW 


DEI 


WRlfE OUTPUT 
WHITE OUTPUT 
WRITE OUTPUT 
DEL ( 1 ) «0«0 
ZNEW ( 1 >°Z (?) 

DO I3ii I»f»t J 
ZNEw())»Z(?«I>* u FL (I) 

WRITE OUTPUT TnR’t Atl32 
WRITE OUTPUT . 

WRITE OUTPUT lARE At 124» (H (2*1 ) tZNEW II) » 1*1 ♦ J) 

FORMAT (2&H0 
FORMAT ( 109H 
1 P 

1?2 R'Ohmb i ( 22RI0 
1?a Format (3(Ei6.t>»fci3.sn 
1?3 FORMAT (1U9H H 

) R DEL 

DFLaVERTlCAL UtR LtCT ION 

U ( I ) aPAOlAL displacement 

SRDbDTrECT SThtbS IN RADIAL DIRECTION 
SCOaOTWECT STRtSS In ClHCDMFtHENTT aL DIRECTION 
SRB«BEmDInR SlKtSS IN radial direction 

SCRaRE'iOlNO SlxtBS IN cIHLUmR ERENTI aL DIRECTION 

HlNTaH 

H«T 

DO 150 T»?tN 

C2«A7 U?* (P-P ' * (R ( I > **2-RU**2) / <2.*CVP) ) aPHR ( I > *UZ ( I ) 

U<I)«A*CVP*(0uU>-C1-CH)/m5(1> 

CI*A5(t) # R(H^LVP 

SRCM j > ate tl ) ♦uZ ( I ) *cVP* (P-PI* (R < I > **2-R0**2) / < 2 t*CVP> ) ) ZC2 
SCO (I » • (DO ( n -Pl*R ( I > *02 < I ) > / IH*Cll 

H2aH*M 

RMOMT ( I ) a-D* ("R U)»XNU*R (1)/H(I))/CI 
CMOMT ( I ) a«D* (R ( A ) /P ( I ) ♦ANU*OR ( I > ) /Cl 
SRB(I)»At*RROf 1 (11/H2 
SCB(I)ab#*CMO- I ( I ) /M2 
FM(naO(l)/R(U 
HO CONTINUE 

1 ARt htlHO 
|«Rt AtlHl 

1 «Rt b<103* U*RU) tU(l> tSHD(T) tSCU(I) tSRMU) tSCH(I) ♦ 


WW 1 TF. OUTPUT 
WRITE OUTPUT 
WRITE OUTPUT 
WRITE OUTPUT 

lpt ) 1 FORMAT ( 3 *>Ht) hmU- 1 « 1 . DISPLACE HENT uNU STRESSES) 

1 x 1 FORMAT ( 1 0 1 H LUCaHON RAlitDlSPt 

ITAM.DlRtSTRFSS RAU.HENOtSl rESS TAw.HENOtSTHFSS) 
FORMAT ( 9 BH A *(T> U(T) 


1«2 

1 SCO(I) 

H 3 FORMAT (l 5 tAE'l*'tb) 
WHITE OUTPUT TARE 
WRITE OUTPUT I **Rt 
IN) 

1 b 5 FormAKB/ho A 


RAUtUlRtSTRESS 

SRU(J) 


ShM(I) 


SCH(I) ) 


At IBS* 

At lHb* ( 1 tR ( I ) tR'UJMTU) tCMOMTU>»FH<T>* F(T)tI* 2 * 


A (I) 


RMOMT(I) 


CNUMT ( 1 ) 
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AUG, OF HUT,) 


1 mOH, FORCE 
1 flb FORMAT (I5«6H'>'6> 

HwHJNT 

80 TO (30n*lMn*iPfi) *NCOM 
1<*H 00 TO tl9(l*2Au»JU0) *n©f: 
IRQ COfcl <1 ,KOMT> b F <*> 

COE ) (?*KU'iT)»vl4» 

COE 1 <3*K0MT)«»F to) 

COE) |4«K0NT)»iitNI 


J«N/? 

COE ) ( S • KOMT ) »i U ( J ) 

IF(KOUT-b) 
lo2 K0NT»1*K0NT 
00 TO 86 

1h 5 WHITE OUTPUT l«Pt 6,198 

white output Impe 4*199 

1q« FORMAT (38H0 INFLUfcNCt COEFFICIENTS OF PLATE A) 
lq9 FORMAT <73 h0 Mi hi M? 

i P FI) 

WHITE OUTPUT 1 w Pt 6* <>1U* c (COEt (I*J) * l*fe) *1»1*S) 
610 FORMAT (6E)0.b) 

80 TO (300»S>*U*0) »NCOnT 


8?0 N6F.»2 
PLSwl.O 
PIS«0.0 


80 TO 6 

240 C0E2U *KONT)»f <*> 
C0E2(2,K0NT)*o U) 
COE? O*K0NT) *E tH) 
COE2(A*KONT) a OlU) 


J»N/2 

C0E2(R*K0WT)». FLU) 

IKKOOT-6) 24f*c:45*24b 
242 KONtbI^KONT 
©0 TO *5 

245 WHITE OUTPUT TAPE b*248 
24© WHITE OUTPUT TAPt 5*199 

248 FORMAT (38H0 INFLUENCE COEFFICIENTS OF PLATE R) 
WRITE OUTPUT TAPE 6* bib** I <C0E2 ( I * J) »<J«1 *6-) * lol *5) 
80 TO (300*300*249) *NCONT 

249 00 250 I«l*4 
DO 250 J»1*A 

250 COE ( T « J) sCOEl U *J) ♦ C(>£2 ( l « J) 

00 255 J»1*A 

PH ( J* 1 ) ®COE I (u*P) ♦COE2 (0*5) 

256 PH(J*?)bCOE1<J*o) + COE2(J*«>) 


0S»1 *0 
J«4 

M*XSImE«F(10* J*2*CoE*PP*OS>*EH) 
00 TO (268*261 *cb2) *M 

261 white nUTPUT TAPE 6*263 
80 TO 267 

262 WHITE OUTPUT TmPE 6*265 
80 TO 267 

257 STOP 555 

243 FORMAT ( 17H0 U*tk/UNDERFLOW) 
245 FORMAT (16N0 CUt is SINGULAR) 
248 BCla-COEtl *1) 

BC2»-C0E(2*1< 

BC3»»C0E (3,1 ) 

8C4»«C0E (4 * 1 > 

NCONT«I 

N0E-4 

H«HINT 

OU 830 I»l*4 
830 80(I)«-C0Etl»2) 


h? 
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60 TO h 

3n0 60 TO (30PtfaO^*bUl»,30))*NUF 
3nl N»«P 
NflEal 
PlSa-PLS 
P JL«-PTL 
K»N/2 

OtU?»OPI.«K) 

NHL»1 
60 TO 838 

3ni 60 TO (808*808) «Nk 
Bn 2 KMAxbM /2 

0EL1-0PUKMAX) 

NRL*2 

PLS»-PLS 

PIU«-P!L 

N0€a2 

8C1»-HC1 
HC2a-HC2 
BC3»-PC3 
BC4»-8C4 
60 TO 83S 
&nb kmAx»n/2 

0fcL2*DFU (K M AX) 

NHL»3 

NR»2 

60 TO 835 
4p5 KMAX»N/2 

OEUlaOFL (KMAX) 

NR»1 

NCONT-3 

NRL»4 

8*5 KmAx«M/2 

00 83ft k«j «km<*a 
636 DSP(K.MR*lK)«uM. (M 
00 838 I*2*N 

JS(I*MR*UK)eOU) 

ROT(I.NH*LK)«P U> 

R0S(t*NR*LK>»S80U) 

C"'S ( I * I 'IR»LK) »!>'*0 ( 1 ) 

RBS(I«4R*t_K) *£>*»< 1> 

CBS(I*MR*LK)«80tim 

RMMT (I*NR.LK)«RmOMT (I) 

CMMT ( I *NR*LK) atrtOwT (1) 

838 FHS(I. <lR*LK)»TOtl) 

60 TO (e>»6*5Qb*bOt>) *NRL 
5ob 0EL»DEL2-0EL1 
SK»ABSP(1./0EL> 

60 TO (840*882) *U 
B 4 O SKPsSK 

60 TO 84S 
842 SKPaSK 

845 WRITE OUTPUT TAPE M410*UiLl 
WRITE OUTPUT TAPE 8*411 *0EL2 
WRITE OUTPUT tnPfc 8*412 *OtL 
WRITE OUTPUT TaHE b*413*S8 
4 ) 0 FORMAT (20H OEM.. OP PLATE A "Flf. 
41 1 FORMAT (20H OtEL. OF PLmTE B "ElB. 
4)2 FORMAT (20H Ott'L. PER CUN** *F15* 
4)3 Format (20h spring hate apis. 

844 60 TO (850,870) »LP 
880 LK»2 

8*0 60 TO (86?*8b2*Pb5) «NCOl*T 
8A2 BC1»-PC1 
BC2«-HC2 
8C3»-HC3 


122 


ut as tr u* 


HG4*«RC4 
PLSwO.O 
PIU«->.00 
NttEwl 
MWNJNT 
NR»1 
00 TO R 
8f& eeiiMixn 
BC2w8D<2> 

BC3wHU<3) 

BC4wBD|4> 

PLSwO.O 
PIl-1.00 
NCONT» 1 
NGE»4 
H«HINT 
00 TO 4 

870 PLS«-<SKF/8KP)*RPlO-0L*CPt*8H*!?KF 
KMA )( iN / 2 

WHITE lUTPl'T IMHK m872*PRL0*PL*^LS*SKF*SKP 
872 FORMAT •'•PRESSURE LOAD ■£!?.&♦ t 

|32H PrFSEi F HtF ItNWT't CHANGE »tl2,S, f 

?32H REACTJVN AXIAL LOAD *E12»5» t 

340H SPPINu LUn$TANT FUR AXIAL LOAD *E1?.5* / 

440H SPrIN(> CONSTANT Fur PRESSURE LOAD •E12.S* t > 

DO BVB NR«»1*2 
SN»1#0 
PL»SM«PLS 
PJ«PRUn»SN 
00 TO (8TB»b?R» *nR 
870 WRITE OUTPUT t«Pfc B»B7R 
00 TO 8«0 

870 WRITE OUTPUT >*Rt R»0^ 

87b FORMAT (//lbH *KOR PLAtE A ) 

87V FORMAT (//lb* »FOH PLA1E « ) 

8P0 00 BB? Kw?»NM*< A 

8P2 DFL (K) wPL^OSPIMNH * 1 > ♦Pl*USPlK»NR»2) 

00 8V0 I®2«N 

U ( I) »US U»NR»i > *PL*US ( I »NR«?) *PI 
SRO ( I ) aROS ( I ».iK» 1 ) #PU*RUS 1 1 » NR , ? ) *P I 
SCO ( I ) aCOS ( l 1 ) APLtCUS ( I»NPt 2) *PI 
SRH ( I ) aR»S ( I *’<R * 1 ) *PL*WuS ( 1 , NR, 2 J *P I 
SCR ( I ) aCBS U > *PL*CUSt I«NR,2) # P1 

RMOMT ( I ) aPMMT t 4* NR, 1 ) *PL7RMMT (| ,NH*2) *PI 
CNOMT <•! ) »CPR t ( I »nR, l ) *PU*CMmT (i,NR» 2>*PI 
Ph ( 1 ) »FHS ( I * N** » i ) *RL HS ( 1 ♦ nR ♦ ? ) *P l 
8Q0 F(I)*R')T(T»NR»1)«PL*R0T (I»NRi?)*PI 
SMAXaO « 0 
UU R9? I»?,N 
DO HV? L B 1 »2 
SN«g*L-3 

SI»SR0(I>7SRHU)*SN 

IF (SMAX-ARSF > /01,^.02*702 
?a 1 SMAX«ARSF IXl ) 

I M A X ■ I 

7r2 S2aSC0tI)*SCHU)*Sli 
S3«ARXF<S1-S2> 

ToJ IF (s^AX-ARbF u*:) ) 70w,fOS,70S 
SMAXaa^SF IS2> 
jNAXat 

7"0 IF ts‘iAi>..«*3) 

UrI S M AX»s< 
l*AXa t 

092 ConT!*''IF 

WHITE HlIPI'T ( «RL It » H7 i * S'' "AX » R ( I MuX ) 

0'JJ FORMAT i-Ax1«oM ST H F-SS tUTENslT^a 


4 
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,;VM AT l 0 * •' 1 ’ 

»M!f "'IHIM -'.M hr 
WlOTF ‘KIHMII ^ ¥' 

i*h n r oiifi'Mi • ,ir 1 •' * J f ** 

wwttf ouIont * ,,e * rv 
wwttK 011101.1 


; ,< , i (.n • uh ** f 
I «t t- ((if? 

Hn <" «P*n 


SiS S 

IN) 

SO 11 070 
lOnU CaU. K*M 
fcND 


SUHHOOTINF 8E,,t V^i? J Y M Jiu , lM , pAT»rX FO»K 

CAUCUl ATION OA INHUT UAU 
SUH A OUT INF SEi.UoL 

UinFNstON 

TANM<)»5.tF'F (»»/COSF(*> 

5 ^^ 11 »sinfuo>-c( 2 ^)*cosF(ao)-c<z»h> 

00 10S i«*id 

. C(id>Bno-c<l*d ... modus 

SilJlifflKiiiiS*"'-' 

GO 10 10V 
S 00 10H I 8 ?d , 
a cti* 2 »*z H +cd*F> 

<* sp»po«pi 

ClU2l a0,n . „ . 

ASGbSP** 2 -P t ** 6 

AGGsSO’TF < ft60) 

AOGbPT/AGG 

AGbATANFIAGG’ 

go 10 (U 0 *U"* i3u> ♦ NQt 

) U 00 11?) I°1* 7 

KOLOaClId) 

20LD»C(I*?> 

21*7010 

’ollrnSF UQ) 4 ? 1 *SINS <A<M 

S! :i!3iIcS5 u5Ui*»*** •*” 

li5 CONTI?" IE . , . 

A6L1 *AT ANF (C 1 1 *** ) ) 

■ AOLTbAtANF (t( M**) ) 

CtUA)BlANFlAoUi-uO 

Ct f«4)»‘TANF<ft"Uf-Ari) 

GO 10 (Un.UdiJOt.NOE 

11 1 WRITE OUTPUT MllH H 

IS FORMAT (1VH0 PtU.OW 

GO 10 1?3 ... 

ipU WRITE OUTPUT S’ 121 

1P1 FOHMATUVHO PtLI.OW 
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i 




i 


AlmAli. lK (ft ( 1 »*• I I 

c ( i *•>> »Hi<r u ••'••tusfr iad) 

DO l?'i )*?*)> 

A<*7bATAN 7(C l I* *1 I 
C(h*'i»»wW«»f. < (A w '* 

C ( 7 • ^ ) *WO*r </«-’•> 

C(l*3)oOtO 

C(7*3)«<i.n 

WWlTF Olllpin I*)*!* Mirth 

UmtF SIRS '•'*■1. **12«* 1 ‘Cl l*d> *1*1*7) 

filiMl? I. HU1 oitMA DF *F*Mt HEtLM) 

FowDAf (Vh ht»*«)AH 2C**l»H 

\ ► i ** ibh wi«) 

lptf FOP* AT (BFlB* »' 

HO SEtUtM 

END 


Iph 
1 Pi 


AI**1BH 


cSImion w'TIncT.ON MU. .FI .NO Mltrtl „ 

SUDWOUTtNF YH mMF tG*N»H»ANU*A»0»NHCl*NBC2*SV»P» l» 

DIMFNS!0N HC Y(ntll('8>» »UO*l« KUM)* OZ(IBO)* OOZC \ M [ ! .! i J*? 1 * 
DIMFNMON 00 ( 1&0) ♦ A1U9Q)* AgU&O). A3 { & BO > * 

A4U*0)« Ab(ibO>* AMIS»0)* A7(150> 


F« 0* UF» RG* 


CyP* Ail A2* A3* AA, AB* A6. A7* 


COMMON* K» t'2* ut‘2* 

l CAUUUtlON| f OK t A- JO H e*C£PT B*Ct 

J*2*l^** 

00 A07 1*1*0 

00 A07 K*).J 
An 7 Y(1*K)«0.0 

y i 9. t ) * I U .«3»*H*Al (2) 

Y (2*2) *“lB»"l' i »**)'* A 1 (2 ) “1*! • am?* A 2 (2) 

Y (2*3)*-A,*lB.*f'*Al (2* 

Y(2*A)*lA.-b.*f'*Mt?) 

Y (2«B) »“b* ♦ (■**»« i <2) 

Y (2*b) *1 • .... 

A1& Y(N»N-4)*1 , , 

YJNoN-^n-bo-’ 1 **! < fl ) 

Y(N*N-2)*lA.»b* # H*Al ( N » 

Y ) o-A »"1 0**H*A1 ( D ) 

Y(N,N) .-IB, ♦1".*M»A\ (K)-12.# m 2*A2(N) 
Y(N*N«*))»)0.M.*H*a1 (N) 

Y (N*2*'*i*l ) •*l* » # * , rt*AA (M /O 
8(N)*-t2.*H2*nS»t^)AP/l> 

A?2 YtN*3*'l*2>*Y(2*l| •u»AAi(9l 

Y(N*3*N*ji*“n*"io »*H*Al 12) •12»*H2 A3 *2* 

Y(N*3»M*A)*Y(m3) 

Y(N*3«D*t>)*rU »A) 

Y (M*3* *1*8) *Y (<•*>) 

Y (N*3oM?> ■) • 

] *PI * (02 (?) *<■' (*) »on?<2) ♦#! (2)#«<2)*02(2) > ) 

A9A Y (?*M*1 «2*N“3l *1 • 
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Y (?«N* \ »***' • 

; j ”3:; 1.. *-1 

y |?#M* | * r* ♦ r’ 1 »» I' »l *U 

rt | ) (no I'll *p*c I *("7 Oil ♦»* (u I (M *wl d>' »*!■•/ • ' 

Oft ft'jfi 

K»I*N* 1 

Y ( 1 1 1»JM bh«A I II I «l « 

V( < I [I 

Y(l«|)»-Jo.«l« «*Mf »«<* ( 1 1 

Y< l*l*Y • a,|W l 1 (1 1 

Y ( I * I ♦<> I ■»l • I HI 

Y < I*K> I ) /I* 

Y ( K ♦«»!> I • Y < 1* i"< I 

Y I K |K«I I "Y ( I * I “ 1 1 

YtK»K*ll*Y(l*l* 1 » 

Y(K«K*?)"V(l*i^*l 

\ ♦Pi*(wn>^ui*i.nz<i>*»ia>*H(i>*u^n»i 

WHITE TUTP-.M IMI'E B^IO 

... rjas^as ’rt.taiWuunj ««" *- zi:\ 

21 1 EoH^ ft T (3PH0 u MatpI* UUCUWATHJ tXCfcpT r OP «#C»l 
00 TO *‘*eH«<*62) •NMC1 

Y (1*2) ■-lO.*l»' • *MMU*M/»12> 

Yll(3l»lB. 

w ( 1 *41 »-b. 


> 11 ) 


**& 


4*B 


♦76 


0 ( I |b!*!N*H*A‘> U) *P (2) *bCi/ (U*CVP) 

00 TO (4t> t i*4br. *4b«.*4?!j) lN tt Cl 

Y ( N*2 * l J*3) ol • 

b(N»2)aHt?)*Bt:<: 

80 TO 403 
Ylli2lil. 

0(1) «nc 1 

00 TO (4b?*4ft>*4bb*4B2) »Nt>Cl 
Y (N42«N*2) *"3« 

Y(N*2f Y*31*"l* *H*XNU/R (2) 

Y(N72»M*4)b10» 

Y (NT2*NTb) »-b, 

b(N*2)»12l*H* |«b|tl (0 vP*A ) 

00 TO 483 

STOP 111 fvPC tHHOH 

00 TO (4b4*49ij*4*Hi*4B4) *nBC 2 
Y(N*l*N-3)»-l» 

Y (KU 1 ♦ * 

Y (NT 1 *N“1 ) B"1 '1 • 

Y (NTl »Nl*10» 4 l2» # h4XNU/K I'tI 

M (Nt!)B-l2!*M*4b (I.) *H(M> *0C3^ (n # CVP> 
00 TO (407*4H ft *4b2*<*V7l *NOC? 

4At Y (2 *Nt?* 2*M*1 > * l ♦ 

0(2*NtP)»H(N)*B 04 
00 TO S04 
4<jO Y(NtI.m)b), 


442 

4«3 

4«4 


no o*o 


HIM* i ) s*HC 3 

GO TO ^ ffs^Cp 

*0 1 V (£*N*? f rffrMrtf* » » r i • 

V ( • #*#M— 1 ) B,J t 
y i p#N*r*i tf*N) 

V I * J (' « 1 ? # #1*1# Kf || )/ R (Ml 

rt <?»*N*«1 ")?»*"* ••'«/* 7 CO ♦ (Ml «(-»C4 > (CVH*A ) ) 
i ♦Pl#»(N)<*nz i 


*> I d 1 J 
MATRIX) 


(iO TO <|04 

h(\k wHTTr ouTnoT i#*^K 

1*1 "uMMAT C3PH.. MWU.H.AHV ^IMUONfi COrtMUFUO) 

UK M INI * 1 
jw 2 *N*A 

WhjTfi TAPF .1*' 

HF WJ MU ,'j 

mM t T FI <1010111 li'M- 
8)3 FORMAT ( J 4 HII I' 

2l* FUHMAT I t)F 1 * • "* i , ., 

WhTTF OUTPHl T«*L Ii»(Im1*1*J) 

Sll CALL SOLV (Y*i *'«*K *(»*SVl 
CALI ORH IHiNiMI'H 
CALL OFW (MtF»*«U'i.) 
wMJtF MlTPltT iwt't **215 

wkITF OUtFiil l*M At2l 7* 1 1 *U U ) «F ( H *b(l) »I>F U • *Ott (t ) ♦ J*2tN) 
21 b FORMAT (A 3 HC 3 L«LotIO | ' < THt F«U«TlOW, F(l) AND 0 ( I ) ) 

2lK FOHMAT (7®H i m *“ 

1 DF f > 0 > 

21? FORMAT (U*btU.b> 

HLTllRH 


End 


210 


2 ib 

3*0 


SUHBOUTJNF ZF* Im*'«*M*P« 2 ) 

CALCULATION OK K 1 1 ) ♦ 2 ( 11 * l’ 2 tl)« DU 2 ( 1 ) * 

SUHpOoTlOF ZF' (L*im»M«W* 2 ) 

StNSE smlTC* J ueClwCOLAR mwc 
ubsealpl mfllow 

OlMFNStON Ctl«*lU)* Ml5u)* 02(150)* 007(1-0) 

common h » n 7 * muz 
CALCdl. AT 10*' Or Mi) 

n:«*c(i*n 

HiloC ( 7*1) 

2 J»C M * 2 ) 

20 *C ( 7 * 2 ) 

X'VBN-? 

SPBpO-Rt 
HbSP/X I 

sJbN + 1 

DO ?ll) 1«1 ♦ vJ 
* 1*1 

H t t ) ®M t ♦ ( X t **2 • I *H 

CALCdl iTlOM OK A ( 1 ) * 0 / ( i) * 0 L '2 ( 1 ) 

Kb 1 

IK (SENSE SWlTu' 3 ) Fib**!/ 

WHITE 11 ) 101)1 I “Ft **300 

FOWL- AT ( 2 NM) LlHGol AN AHL Htl LON NLAtF) 

DO FI* 1“1 *<* 

C ( 1 • ! ) »C ( 2 * 1 ) 

C (? * 5 ) bC ( ?* 5 ) ♦ •>«**' 
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2l7 DO ?<►< 1 i»l*U 

IF (K-71 214*ZCb*22A , 

2i 4 CVl*AH9F (c (M J> **<;.» (H < U"C (*• 1) ) ## 2« > 

C72*SORTF (CV1 ) 

Z(I)«C(K»2)*C(A*A)#CV? 

80 TO 218 

2?6 Z ( I ) »C (K*?l *C t ^ > * (H 1 1 * CK* 1 ) ) 

2 1 B IF (C (K *S) -R ( I ♦ A > ) ?38*2J6»2aO 
2*g IF |CIK»l* 5 l*»U # in 239*239*237 
2? 7 K*K*1 

80 TO 290 

239 K»K*2 

240 CONTINUE 
J»N-1 

DO 245 I»3*J 

07 (I ) ■ (2 ( T-2) “0**Z ( I"1 ) *8«*z 1 I* 1 ) m l < ) / ( 12. *H) 

DOz 2 2)»ao.*Z(i>-l‘5.*Z(‘i)- < *.*Z(3) ♦H,«Z(4)-6.*Z (8) ♦/ (fc{ >/ 

DZ <n)«<-Z(N- 3) ♦b»*7<N-2>«18**7 (N-l > *1 0,*Z (N) *3.*Z (N*l ) > f ( 12* M> 
ODZ (N) * (Z (N*>9 ) “0»*7 (N"3 ^ 4 i A«*7 (N-2)-9.*Z (N-l ) -1 5»*7 (N) 
1*10,*7<NM> )/U4.*H?> 

WRITE OUTPUT ImFL N.2bO*H 
WRITE OUTPUT l«Mt M«2bl 

WRITE OUTPUT I^F'k 0.252* ( I tR ( l ) *Z t I ) *02 ( I ) *007 < I ) * I*2*N) 

2A0 FORMAT (35H0 INTERVAL M/E Op BELLOW PLATE H»E13.5> 

2nl FORMAT (4H i*iU*i H<l)*15H 2 ( I ) ♦ 1«>H 

11.17H Ul'7 ( 1 1 ) 

2E.2 FORMAT (Ib»4Flb»b) 

RtTUR" 1 
EMO 


DZ ( I 


C 


CALCULATION! Ok LOt F I C lE'^T Al 
SURR0UT1NF AF 1 ' lH./*N*H* Xi»U* U « A ) 


OIMFNrION h( 13U)» 7 ( 1 *>o ) • DZ ( 1 50) * ODZU50)* AJ <150) • A 2<1^0-* 

I 6J(ibu)« A9(i5u). AMl^O)* A6(l50)« A7(150)* F(150)# 

. 0(130) * t)P(lSO). UOUbU) ' 

COMMON k« Dz* RD/ * F. R. OF* 06. CV**# Al* A2t A3* A4. *5. A6» A7* 


) ANU* ■!« A 

J«N*1 

Of) 305 1*1 .J 
CVR«0. 159155 

ZPl*l.>0Z(H**2. 

ZR2*0/ t I ) *r>DZ UJ/ZPl 


RC* 1 • /'* ( I ) 

ZR3=S()rTF (1«* •Z(l)*«2») 

Al (1)*RC-7P? 

A?( 1 1 aRC* (PC*A iD*/PF) 

A3(I>*RC*(RC-» -Li**!-’?) 

A4 ( 1 ) ■JC*07 ( 1 1 

Ab<I>»CVR*R(; # /' fJJ 

ZR4B0/ ( I ) WPC* * 1 «R*l-U/ ( 1) ZZRl 

AM ( 1 loCVRAROW/p'A 
AMl)«CVR*n/ ( I l*RL*XNU 
3*5 CONTt'i iE 

KKfTF lUtPUT ! U k m. 315 
WRITF OUTPUT l<‘ft m.3 1M 

wRIlE HJlRnt i«*'L f> . 31 7 ♦ I f « A 1 ( T ) « AZllI* 
1 ) * A 7 ( 1 ) « 1*1* •'> 

Jib FORMAT (IMHO **"<CllO> Mill 
3lh F i.iWf-’At (bH 1*1' f Al(I)«l3R 


A 3 ( I ) ♦ A4(I) . »5(I) • A(i ( 1 


A? f I ) « 1 3 m A3 ( 1 ) * 13 


oo 


AT' ( I ) » 1 -1h 


AhU) »> 3 H 


AM!) ) 


1 H 

3l 7 FORMAT 
R) tuH'i 
E i,0 


AA (!) *1 Jh 
(Id* h iitb) 


SUHwnuTIw t.'E*' 

SURROMTINF lit.** IM»I'»Z«|U) 

dimension /davit r<Z 1 1 t»o » 

0*<a)»<-3#ViN-l<%«/ <*>♦**•*/ 1*)***** (4)*Z(S) !/(»••« 

0/ (n) * (M"3) *&•* / (N*2> (N»J )♦ i0* # ? (NJ (N^l) ' 

J*>M 

00 70m t*3*J 

7n» 0/ ( J ) * < 7 < I «i!)»'a.*3(|-l)*M.*7(!*l)-<'(I*2))/()?.*H) 

RETURN 

END 

SUHROUT INf IN I CO 
SUHHOilTINF IN I tl) (F*H*D*S) 

UlMtNStON fUSVlt SllSvl 

COHV»0«0 

j»n/2-1 

on ?o *®l*J 

1»2*K 

S(K)*C1NV*H/3. 

C0NUbC 0NV*A.*I‘ ll*ll*FU*il*MI) 

?U CONTINUE 
J8N/2 

S(J)«C0NV*H/3. 

RETURN 

END 


2) A 

2l« 


2 no 

2*2 


3n0 
3) 0 

312 

313 
3l A 

316 

6n6 


SUHPOi'TINF 5 >OlV 

STOP AAA S 1 N >OLAH MATRIX 

SUHPOUTINF 801 » (T«R*N*F *6,D> 

FORMAT (BFIA.'M 

FORMAT 1 19W 0 nftTwi a CHECK) 

DlMF'NSTON Ytio2*l(i2)* KllO?). FC1S0J* 0(150)* F" (150) 

CONST«n.lOOO 

FORMAT (5E) 3*6) 

wPtTE OUTPUT T**K m *202 *CONS f 

FORMAT (33H0 8CAtt FACTOR U8*U) IN YBMATRIX«E13.B) 

CaLCUUTIOM OF Hfck Y* H MATRIX 

J»2*N*2 

00 300 I ® 1 * U 

Bit) «r,0NST*B ( 1 1 

UO 300 K®1,J 

Y(t*K)®CONST*Y U*K) 

WRITE OUTPUT lAFt h,310 
FORMAT (15H0 ' t* H MATRIX) 

WRITE OUTPUT 1AFE 6 * 31 2* <P < I > * I®1 *J> 

FORMAT (8tlA»b) 

IFtsEMRE SWITCH 1) 313, 606 

WRITE OUTPUT ImRE 6*31a 

FORMAT t lbHO ' tw Y MATRIX) 

WRITE )UT PUT TmF'E 6*316* (! Y (I »K) »K*I » J) »I«1 * J) 
F0RMAT(9E13*b) 

M»XSIMEQF ( 102*0* 1* Y*B*D*E> 

WRITE OUTPUT ImME 6 *600 *U 


2,*m) 


<* 



00 TO (6Uri»t>Ul »M 

fenj WHITE HlTPMT 1 
h<) TO Ml? 

t,n£ WHITE OUlPUT l.'Hf, h,hhh 
GO TO ol) / 

MO format |6HA 

bn 3 FOWMATU7HO > vtR/UUOt HH-Uw) 
tin# F OPMAT I lfiHfl i iS StNOljLAh) 
<jn 7 STOP AAA 
MB JaN+1 

DO Fill I»1*J 
Ka I *N* 1 
FIllsTIUI) 

61 l GU)aY(Ml) 

J»2#N*? 

DO 61ft I«1«J 
H ( I ) art 1 0 
E(!)«0.0 

616 

DO 61V I»1*U 
DO 61V K«1 » J 
619 Y(1«J)*0.0 

READ TAPE. 3*Y 
REWIND 3 
DO 62# 1*1 ♦ J 


DO 62# K»1 »J 

6?# B <1 1 »M < 1) aY ( 1 *E. <K) 

WHITE OUTPUT TAPE 6«218 ...... 

WHITE OUTPUT t*Ht ft *21# t ( t> (I ) 1 1«1 * J) 


RETURN 


END 
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APPENDIX B 


DESCRIPTION OF COMPUTER PROGRAM "NONLIN" 


Program NON LIN was developed by Battelle Memorial Institute , 
under Air Force contract. It is generally capable of predicting the 
nonlinear behavior of bellows plates under a wide range of conditions. 

It was used in this study for all nonlinear range analysis and for the 
linear range analysis of the S -shaped bellows plates. 

Because NONLIN is describ d in detail for the interested reader 
in Reference (6) which is available i m the DDC collection, omy a 
brief description will be presented here. 


Capabilities 

NONLIN can treat any bellows plate or shell which can be 
divided into segments having the following basic chapes: 

(a) cylinder 

(b) spheroid 

(c) ellipsoid 

(d) parabaloid 

(e) conical section 

(f) toroid 

The thicknesses of these elements can vary along the meridional 
direction. Thus, it can be seen that virtually any axi symmetric 
plate shape can be treated. 

The possible applied loadings are as follows: 

(a) pressure acting normal to the plate surface (lb /in ) 

(b) weight of the shell itself (lb/in ) 

(c) dead weight acting on the shell (lb /in ) 

(d) boundary forces, moments, or displacements acting 

at the initial and final edges of the structure. 
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Load (a) can bo specified to vary along the radial direction but must be 
axisymmetric unless it can bo expressed by a Fourier harmonic 
series with eight or less terms. Loads (b) and (c) must be constant 
over the element and the shell centerline must bo vertical. 

The elastic properties of the shell (i. e, Young 1 s modulus and 
Poisson' s ratio) can vary along the meridional direction. 


Computational Method 

NONLIN uses a method developed by Kalnins (7, 8). The 
general nonlinear shell equations are reduced to a set of six first 
order nonlinear differential equations. These are solved by a multi - 
seement model which combines direct numerical and finite difference 
methods. This method is coupled with a systematic iteration procedure 

to solve the nonlinear equations. 


Program Limitations 


The program was written in FORTRAN IV for the CDC-6400 
computer. It consists of one main program and 15 subroutines and 
is self-contained except for standard library functions. 


The major inherent limitation is that the total number of 
segments which can be treated is 60. For simple plate shapes, this 
is not a serious limitation, but the number of convolutions for a 
complex plate shape becomes limited. This precludes the study of 
the effect of edge conditions several convolutions away from the 
boundary or the simulation of a convolution in the center of a long 
bellows core. 


Occasionally, in runs for this study where near 60 segments 
were used, the iteration procedure required for the nonlinear solution 
did not converge. 


Other than this segment size limitation, the program was 
found to be extremely flexible, and economic for predicting bellows 
plate behavior. 
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APPENDIX C 


COMPARISON BETWEEN "SEALOL" AND "NONLIN" 
FOR LINEAR RANGE ANALYSIS 


Program NONLIN can also be used to obtain a linear solution 
for the stresses and deflections in bellow. Thus, for the benefit 
of a potential future investigator who requires linear range solutions, 
the following user -oriented information is documented. 

For linear solutions, the required computer time for the two 
programs is comparable, NONLIN probably being somewhat faster for 
most shapes. Preparation of input data is complex for both programs 
and requires comparable time. Here NONLIN is probably also more 
convenient due to the documentation available in Reference (6). 

If the desired solution is for a single convolution in the central 
portion of a long bellows core, SEALOL is likely a more accurate 
tool because SEALOL treates a single convolution whose boundary 
conditions are provided by an adjacent entirely similar convolution. 

It, thus, is not capable of predicting the influence of the end fixtures 
on the behavior of the convolution. NONLIN, on the other hand, is 
capable of treating only 3 or 4 convolutions adjacent to the end 
fixture. Thus, it is capable of predicting end effects but cannot 
represent a convolution far away from the end fixture. 

Both programs are capable of treating a variety of basic plate 
shapes. Because of the method used for describing plate shapes, 
however, NONLIN, is more suitable for treating complex convolutions 
consisting of coupled segments of basic shapes or stiffening rings. 

It also is capable of treating variations in plate thickness across the 
span where SEALOL is not. 

Another fundamental shape limitation for SEALOL which does 
not apply for NONLIN is that it should not be used to treat any shape 
where the local plate inclination approaches 90° (e. g. near the edges 
of a full toroidal span) since the form of the equations used has a 
singular point in this case. Thus, it is suitable only for shallow 
plate elements. 

Finally, the programs were each developed for a different 
computer although each is written in FORTRAN. Modifications in 
control and format statements are often required in adapting a program 
written for one machine to another machine. This was found to be 
the case for NONLIN which was originally written for a Control 
Data computer. Modifications would have been required in the program 
to use it on an IBM machine. Thus, for convenience, it was used 
on a CDC-6600. The FORTRAN II used for SEALOL is compatible 
with the IBM 7090 and 7094 systems. 
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Another factor affecting the uso of program NON LIN is the largo 
momory required, approaching 100,000 words. This again is compatible 
with the CDC 6600 system which has a 100K memory, but is too large 
for those systems which only have 32K available without augmentation. 
Thus, a dimension statement modification would be necessary and the 
number of plates treated is limited. 

On balance, it is recommended that the Battelle program, 

NON I -IN , be used for any case where it is suitable, because of the 
available documentation and its inherent flexibility and generality m 
treatment of plate shapes. 
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APPENDIX D 


TABULATION OF COMPUTER RUNS AND PLOTS OF THE 
STRESS DISTRIBUTIONS FOR THE 
SYMMETRICALLY-STACKED TOROIDAL BELLOWS 


Computor Runs 

A tabulation of all computer runs made on Program NONLIN 
to investigate the behavior of the symmetrically -stacked toroidal 
bellows plates is given in Table D-l. This data is plotted and 
discussed in the main body of this report. 

The nomenclature used in Table D-l is defined in Figure D-l. 


Stress Distributions 

Figures D-2, D-3, and D-4 shows the distribution of the four 
principal stress components as predicted for opposed plates having 
included angles of 180°, 120° and 60°, respectively. The loading 
assumed for these plates was a combined 500 psi pressure differential 
and an axial plate deflection of 5 mils. For each plate, a 5 mil 
thickness was used. 

In Table D-l, the maximum combined stress intensity (twice 
the maximum shear stress) is tabulated for the various critical cases. 
This value is obtained by taking the maximum local difference between 
principal stress components. 
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(f >£-0£ ® included angle of the toroid* degrees 
radius of the toroidal arc, inches (2. 54 x 10 
plate thickness, inches (2. 54 x 10 m)^ 

axial plate deflection, mils (2. 54 x 10 m) 


m) 


-5 


P 

n p 

AN p 
A D 

e 

°mm 

^mb 

O' 

max 


radial plate deflection at inner edge, mils (2. 54 x 10 m) 

-5 v 

radial plate deflection at outer edge, mils (2. 54 x 10 m) 

pressure load, psi (6. 89 x 10 3 newtons/m 2 ) 

axial reaction force per unit circumference at the sealing 

*2 

diameter, lbs/in (5.67 x 10 " newton/m) 

the difference between the actual axial force, N , under the 
applied loads and that calculated from a linear solution 
with the applied loads, Ibs/in (5.67 x 10-3 new ton/m) 

change in effective diameter, percent ^ 

meridional maximum stress, Kpsi (6,89 x 10 newton/m ) 
meridional bending Btress, Kpsi (6.89 x 10 newton/m^ 
combined maximum stress, Kpsi (6,89 x 10 newton/m ) 


Figure D-l. Definition of Nomenclature, 
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Figure D-3 Stress Distributions in 










APPENDIX E 


TABULATION OF COMPUTER RUNS AND PLOTS OF 
SA^PI-T stress distributions for the 
nested toroidal bellows 


Computer Runs 

This appendix presents the data obtained from computer runs 
lade using program NON LIN to perform a parametric analysis of the 
nested toroidal bellows. The definition of nomenclature is given in 
Table £1 Most of the variables are also shown graphically in 

ot the" main body. A tabulation of the computer results if 

given in Table E-2. 


Stress Distributions 

In order to illustrate the nature of the stress ^^utions 
wnich occur in the nested toroidal bellows, the case of Run No. 29 
which is close to the final recommended configuration was chosen 
an* 1 example. The general configuration of the four 

f nr the computer run is shown in Figure E-i. xne sneiis nav 
an included angle^f 160° and the dimensions of the shells and the 
stiffener plates are shown, in the figure. 

The four principal stress components are plotted J^^or" en^ 
through E-5 for four successive plates from the fixed b un Y 
The stresses are plotted across the span versus angle, <f> , which is 
^Ened as increLing from the inner to the outer edge as shown in 

Figure E-l. 

Comparing the first shell to the t third < th ell. He. an be 8ee ” tha * 
some variation iccurs near the outer edge which would be expected 
since the first shell is at the fixed boundary at the outer edge. xhe 

S3SS3 

similar plates of succeeding convolutions beyond the second. 


144 


Table E-l 


Definition of Nomenclature 
pressure, psl 

axial deflection, mllll -Inches 
radial deflection, milli -inches 
span, inches 
mean radius, inches 
radius of toroid, inches 
shell thickness, milli -inches 
length of stiffener ring, inches 
thickness of stiffener ring, inches 
included angle of toroid, degi'ees 
angular deflection, milli -radians 
number of convolutions 

axial reaction force per unit circumference, lbs/inch 

radial reaction force per unit circumference, lbs/inch 

reaction moment per unit circumference, in -lbs/inch 

the difference between the actual axial reaction force 
per unit circumference at the sealing diameter under 
the applied loads and that calculated from a linear 
solution under these loads, Ibs/in (5, 67 x 10 
newtons /m) 

combined maximum stress, K psi (6.89 x 10 
newtons /m 2 ) 


Tabulation of Computer Results 



l 

i 


i 

i 
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Table E-2 (Continued) 
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Table E-2 (Continued! 



















































































































































































































































































si (3.445 j: 10 newtons 



Figure E-l. Geometry for Run No. 29 













•e E-3. Stress Distributions for the Second Sh ell from the Fixed End, 

Run No. 29. 





Moridional Bonding Stress 


Tangential Membrane Stress 


.Tangential Bending Stress 

Meridional Membrane Stress 



Angular Coordinate, 0 , Degrees 


Fieure E-4. Stress Distributions for the Third Shell from the Fixed End 


Run No. 29 




Stress, (psi x 10 " 3 ) (newtons/m x 1.45 x 10 



Figure E-5- Stress Distributions for the Fourth Shell from the Fixetl 

End. Run No. 1 % 
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APPENDIX F 


STIFFNESS RELATIONSHIPS FOR THE RADIALLY 
CORRUGATED BELLOWS 


Introduction 

The nonlinear analysis of tho simple opposod toroidal bellows 
showed that the change in effective diameter over the desired rango of 
axial deflect,’ on is beyond the design goal. 

The normal toroidal plate is very stiff circumferentially. Thus, 
as it is deflected, each plate fiber resists moving in or out with 
respect to the plate centerline. This produces radial bending or 
curvature change in order to accommodate the deflection. The 
curvatui e change or angular deflection across the span is plotted in 
Figure 30 of the main body of this report. The plate distortion 
produced in this manner could have the effect of altering the reaction 
forces and moments produced by subsequent pressure loading and, thus, 
change the effective diameter. 

It has been hypothesized that the corrugation of the bellows 
plate around its circumference (corrugation peak lines and valleys 
lying in the radial direction) would mitigate the plate distortion pro- 
duced by axial deflection and, thus, improve linearity or constancy 
of effective diameter. 

Unfortunately, analytical techniques within the scope of this 
study cannot treat a non -isotropic bellows plate of this type. Thus, 
it is recommended that a radially corrugated opposed toroid bellows 
be fabricated and tested for linearity compared with a similar 
uncorrugated one. In order to arrive at a design for this experimental 
bellows which produces significant variations in the relevant plate 
stiffnesses, a simplified analysis was performed to determine the 
important stiffness ratios as functions of corrugation geometry. This 
analysis is summarized in this appendix. 


Geometry 

The geometry assumed for the analysis is shown in Figure F-l, 
where the principal nomenclature is also defined. For simplicity, a 
flat-topped corrugation was used, However, it is felt that the results 
provide a good approximation for the characteristics of any corrugations 
having a half wavelength, \ , and a mean height, h. 
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- included angle of the toroid 

- plate thickness 

- angular coordinate 

- radius of the b > se toroid 

- height of the corrugation 

- half wavelength of the cor rugation 


Figure F-l. Schematic Diagram and Definition of Nomenclature 




Radial Direct Stiffness 


The radial direct stiffness per unit length for an uncorrugated 
plate is simply 

K = A E = Et 
RD o 

since the area per unit length is 


A = t 
o 

For the corrugated case of Figure F-l, the total area per half 
wavelength, X , is, 


A x = t (\ - h + /Th) 


Thus, the radial direct stiffness becomes, 


and the correction factor or ratio of corrugated stiffness to 
uncorrugated stiffness is given by, 


( k rd^ 


RD 


K 


corr, = 2 


+ (/T - i) f 


RD 


Radial Bending Stiffness 

Neglecting Poisson effects as an approximation, the radial 
bending stiffness for the uncorrugated plate is, 


K 


RB 


I E 
o 



Taking the bending moment of inertia of one half wavelength of the 
corrugated plate about its midplane gives, 


(F-l) 

(P-2) 

(P-3) 

(P-4) 

(P-5) 

(F-6) 
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- + (X -h) « (-I-) 2 ♦ (4-) 3 (i> 

“ "TT | ^ " h * + [ 3 ^ " h ^ + ["It] | 


Thus, the correction factor for the corrugated case becomes, 

LA - p ^ i r, n 2 


RB 


1 — * i - £ + 3 [i - (i - 4*) x] [t] 


Circumferential Direct Stiffness 


For the uncorrugated case the circumferential direct stiffness 
is the same as the radial direct stiffness or, 

K = Et 
CD 

For the corrugated case, one must consider the quarter wavelength 
free body diagram of Figure F-2 where symmetry dictates that the 
corrugation midpoint is equivalent to a fixed boundary and that the 
midplane is an inflection point with no moment. At the corner of the 
corrugation segment, equilibrium dictates that, 


M l aF t 


and the deflections can be written at, 


6 = L 1 L AL 

®1 2 Et 

M 1 _ Fh (\-h) 


2 El 


4 El 


(F-7) 


(F-8) 


(F-9) 


(F-10) 


(F-ll) 


1$8 



A 

\\\\\\\\\\ 

T 

h/2 

i 

midplane — - " * " 




V- 


\ 

\ 

\ 

\ 

\ 

\ 


corrugation 

midpoint 



Figure F-2» Free Body Diagram of a Quarter 
Wavelength of a Corrugation Under 
Circumferential Load 
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where 


(F-12) 


X = 


12 


Defining 6 2 as the deflection at the end under F when the corner is 
fixed gives, 


5, ■ -£* (4 -> (t-i + 


1 1 1 


2 -jj 'eF 'X' - 7T< yi' J 3 EI j-j 

The total deflection, 6, can be obtained by adding the effects of the 
above .according to, 


»• \ + e i T +6 z 


I (fx-h h 1 [ h 2 (X -h) h 3 _ 1 J. 
= -2E j[— ^ t J + [ 4 b/T -\ Ij 


Thus, the equivalent direct stiffness is, 

F X 

( K CD^corr. ~ ~ 


= ET 


j 1 - (1 '7r )( ^ ) + 3 1 ' (1 Wj 




and the correction factor is 


f CD = ‘ ^ 


- 7? ^ + 3 [‘ 




2 \ -1 


(F-13) 


(F-14) 


(F-IT) 


(F-16) 
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Circumferential Rending Stiffness 


The corrugations will not have an effect on the circumferential 
bending stiffness. Thus, the correction factor is simp y, 


£ CB 


1 


(P-17) 


Numerical Values 

If it is assumed that the corrugatiomampUmde va^e. acre.. 

the span from aero to a maximum amplitude, H, at the tmap 
according to the relations, 


0 - 0 

h = 2H ( 


4 > 


■) 


9 . 

h = 2H ( " 


<*> 


-) 


0 < 9 < ~ 

o 4 


f< 9 < 9 i 


(F-18) 


*TT 

then the following values can be found at the midpoint, 9 - 2 ’ 

the four correction factors, 


RD 


RB 


CD 


CB 


= 1.083 

= 17. 5 

= 0.056 
= 1 


where 


H 

t 


= 10 


mils (2.54 x 10 “ 4 meters) 


= 50 mils (1.27 x 10" 3 meters) 


= 4 mils (1.016 x 10 " 4 meters) 


unoorrugated eise should be sufficient for experimental evaluation. 
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Use of the various relationships of this appendix also allows 
prediction of the effect of any other corrugation configuration of this type. 
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